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(ABSTRACT) 

Layerwise finite element analyses of geodesically stiffened cylindrical shells are presented. The 
layerwise laminate theory of Reddy (LWTR) is developed and adapted to circular cylindrical 
shells. The Ritz variational method is used to develop an analytical approach for studying the 
buckling of simply supported geodesically stiffened shells with discrete stiffeners. This method 
utilizes a Lagrange multiplier technique to attach the stiffeners to the shell. The development of 
the layerwise shells couples a one-dimensional finite element through the thickness with a Navier 
solution that satisfies the boundary conditions. The buckling results from the Ritz discrete analyt- 
ical method are compared with smeared buckling results and with NASA Testbed finite element 
results. The development of layerwise shell and beam finite elements is presented and these ele- 
ments are used to perform the displacement field, stress, and first-ply failure analyses. The layer- 
wise shell elements are used to model the shell skin and the layerwise beam elements are used to 
model the stiffeners. This arrangement allows the beam stiffeners to be assembled directly into the 
global stiffness matrix. A series of analytical studies are made to compare the response of geodes- 
ically stiffened shells as a function of loading, shell geometry, shell radii, shell laminate thickness, 
stiffener height, and geometric nonlinearity. Comparisons of the structural response of geodesi- 
cally stiffened shells, axial and ring stiffened shells, and unstiffened shells are provided. In addi- 
tion, interlaminar stress results near the stiffener intersection are presented. First-ply failure 
analyses for geodesically stiffened shells utilizing the Tsai-Wu failure criterion are presented for a 
few selected cases. 
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Chapter 1 
Introduction 


1.0 Background 

Laminated composite shell structures have found varied applications in complicated 
aerospace structural systems. This is due primarily to the advantageous properties of 
composite materials such as high strength-to-weight and stifTness-to-weight ratios for 
weight sensitive applications. Additionally, composite structures have a high fatigue life, 
corrosion resistance, low fabrication cost, and are tailorable to the loading environment. 
Aerospace applications using composite structures are almost limitless, but often require 
the use of sophisticated analyses to determine the response behavior to external loads. 
This is because laminated composite materials consist of two or more layers that are 
bonded together to achieve desired structural properties. Material properties of lami- 
nated composites are discontinuous through the thickness because of the different ma- 
terial layers in the laminate. Thus, the analysis of composite structures is quite 
complicated due to material discontinuities across the laminate interfaces, bending- 
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stretching coupling in the laminate, and the geometrically nonlinear effects. Traditional 

analysis methods applied to isotropic materials cannot be applied directly to composite 
materials. 

As new applications of composite structures evolve, so also the analytical techniques to 
study these applications must also evolve. Existing metal aircraft design methods permit 
the skin panels of some structural components to buckle under various loading condi- 
tions. Hence, these structures are designed to have postbuckling strength. Before com- 
posite structural components can be designed with similar buckling response, their 
strength limits and failure characteristics must be well understood [1,2]. Grid-stiffening 
concepts based on new, automated manufacturing methods such as filament winding 
where the co-curing of stiffeners and skin is achieved hold great potential for cost 
savings. Additional applications of stiffened shells may be found in aircraft fuselages, 
rocket motor cases, oil platform supports, grain silos, and submarine hulls. 

Accurate design analysis of stiffened circular cylindrical composite shells is of great im- 
portance in the aerospace industry as it relates to aircraft fuselage design. The objective 
of this study will be to concentrate on the analysis of geodesically stiffened cylindrical 
composite shells subjected to compressive loads. The analysis will include a study of the 
stiffened shell buckling and stress analyses. See Figure 1 for a description of the 
geodesically stiffened shell system. Most previous analyses of stiffened composite shells 
have utilized either a smeared stiffener approach or a linear finite element analysis to 
determine the buckling loads. Although few, nonlinear analyses of stiffened shells are 
typically performed using the finite element method. Analysis of stiffened composite 
shells must include the failure characteristics of the shell structure including general in- 
stability, local stresses, interlaminar stresses, and failure analysis. 
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Traditionally, in order to capture the localized effects in laminated composite shells a 
three-dimensional (3-D) finite element analysis must be used. Further, a fully nonlinear 
3-D finite element analysis must be performed to characterize the structural response in 
the postbuckling regime. Unfortunately, if a laminated composite shell is modeled with 
3-D elements an excessively refined mesh must be used because the individual lamina 
thickness dictates the aspect ratio of the elements. The aspect ratio of the elements must 
be kept reasonable to avoid shear locking. Even in localized high stress regions a 3-D 
analysis will be computationally intensive and expensive. 

The motivation of this research is to develop an accurate analytical methodology for the 
study of stiffened circular laminated composite shells without applying a costly nonlinear 
3-D analysis. The analysis should be accurate in the nonlinear region and provide for 
any localized high stress regions. The interlaminar stresses near the stiffener intersections 
of stiffened structures is of interest to shell design engineers. Moreover, the efTects of 
these interlaminar stresses on the structural integrity of stiffened shells has not been de- 
termined. The literature review in the next section provides a background for this study. 
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1.2 Literature Review 


The purpose of this literature review is to present the current state of analysis of stiffened 
composite cylindrical shells. Also, included are discussions of shell theories, finite ele- 
ment methods, discrete stiffener approaches, and failure mechanisms in composite ma- 
terials. This should provide sufficient background for the detailed theoretical and 
numerical work which follows. 


1.2.1 Shell Theories and Finite Element Applications 

The first classical theory of shells was proposed by Love [3] in 1888. The basic premise 
of Love's paper is the KirchhofT-Love theory in which straight lines normal to the 
undeformed middle surface remain straight, inextensible, and normal to the deformed 
middle surface. As a result, the transverse normal strains are assumed to be zero and 
the transverse shear deformations are neglected. Love's theory can be applied to thin 
shells where the shell thickness is small compared to the least radius of curvature. An 
improvement to Love's work was made by Sanders (4) when he presented a theory to 
remove the strains for small rigid body rotations which are erroneously predicted by 
Love's theory. The thin shell approximations of Love requires that the thickness of the 
shell is small compared with the nominal radius of curvature. Donnell [5] removed the 
thin shell approximation of Love by developing a theory for shallow shells. Reissner [6] 
and Mindlin (7) each developed shear deformation theories for plates and Reissner 
[8,9,10] extended the theory to include transverse shear deformation in shells. Surveys 
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of classical linear elastic shell theories can be found in the works of Naghdi [11], Bert 
[12), Krauss [13], and Fliigge [14]. 

The use of composite shell structures has forced the development of appropriate shell 
theories that can accurately account for the effects of bending-stretching coupling, shear 
deformations, and transverse normal strains. Ambartsumyan [15,16] developed the first 
analysis that incorporated bending-stretching coupling. Ambartsumyan's work dealt 
with orthotropic shells rather than anisotropic shells. Dong et al. [17] developed a theory 
for thin laminated anisotropic shells by applying Donnell type equations to Reissner's 
and Stavsky's [18] work for plates. Flugge's shell theory [14] was used by Cheng and 
Ho [19] in their buckling analysis of laminated anisotropic cylindrical shells. A first ap- 
proximation theory for the unsymmetric deformation of nonhomogeneous, anisotropic, 
elastic cylindrical shells was derived by Widera et al. [20-22] by means of asymptotic in- 
tegration of the three-dimensional elasticity equations. The laminated shell theories 
discussed thus far are based on the Kirchhoff-Love assumptions and therefore are only 
applicable to thin shells with mild material anisotropy. Application of such theories to 
layered anisotropic laminated composite shells could lead to as much as 30% or more 
errors in deflections, stresses, and frequencies according to Reddy [23]. 

The effects of transverse shear deformation in composite shells were introduced by 
Gulati and Essenburg [24], Hsu and Wang [25], Zukas and Vinson [26], and Dong and 
Tso [27]. The development presented in [24] is based upon the shell theory given by 
Naghdi [28,29] and assumes symmetry of the elastic properties with respect to the middle 
surface of the shell. The theory presented in [26] also includes the effects of transverse 
isotropy and thermal expansion through the shell thickness. The theories of references 
[25,27] are only applicable to orthotropic cylinders. Whitney and Sun [30,31] developed 
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a higher-order theory for laminated anisotropic cylindrical shells. The theory includes 
both transverse shear deformations and transverse normal strain. Reddy [23,32] extended 
Sanders theory for doubly curved shells to a shear deformation theory of laminated 
shells. The theory accounts for transverse shear strains and rotation about the normal 
to the shell midsurface. Reddy and Liu [33] proposed a higher-order shear deformation 
theory for laminated shells. The theory is based on a displacement field in which the 
displacements of the middle surface are expanded as cubic functions of the thickness 
coordinate, and the transverse displacement is assumed constant through the thickness. 
This displacement field leads to a parabolic distribution of the transverse shear stresses 
and therefore no shear correction factors are used. Librescu [34,35] developed a refined 
geometrically nonlinear theory of anisotropic symmetrically laminated composite shal- 
low shells by incorporating transverse shear deformation and transverse normal stress 
effects. The theory was derived using a Lagrangian formulation in which the three- 
dimensional strain displacement relations were modified to include the nonlinear terms. 

Recently, Reddy [36] developed a layerwise laminate theory which yields a layerwise 
smooth representation of displacements through the thickness. The layerwise laminate 
theory of Reddy (LWTR) reduces the 3-D elasticity theory to a quasi 3-D laminate 
theory by assuming an approximation of displacements through the thickness. Reddy 
[37] and Reddy and Barbero [38] extended the LWTR to the vibration of laminated cy- 
lindrical shells. Further study of laminated shell theories may be found in papers by Bert 
and Francis [39] and Kapania [40]. 

A large number of different finite elements have been formulated for the static and dy- 
namic analysis of isotropic and anisotropic shells. One of three approaches are usually 
followed in shell finite element theoretical development. The first approach involves the 
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development of finite elements from an existing 2-D shell theory [41,42]. In the second 
approach, 3-D elements based on three-dimensional elasticity theory are used [43,44], 
For the third method, 3-D degenerated elements are derived from the 3-D elasticity 
theory of shells [43-49]. One of the earliest uses of finite elements in layered composite 
shells was provided by Dong [50] on the analysis of statically loaded orthotropic shells 
of revolution. Other authors [51-54] continued the development of finite elements appli- 
cable to laminated composite shells. Nonlinear analysis is critical in the study of shell 
structures. The nonlinear response of shells under external loads was published in refer- 
ences [31, 47-49, 55-59] among others for laminated composite shells. A more detailed 
discussion of laminated shell finite elements may be found in [40]. 


1.2.2 Structural Analysis of Stiffened Shells 


The circular cylindrical shell is used extensively as a primary load carrying structure in 
many applications and is therefore subjected to various loadings. Design limit loads of- 
ten result from general or local instability due to the action or interaction of pressure, 
axial, torsional, and thermal loads. The elastic stability of monocoque isotropic cylinders 
is well documented in the open literature [5,14,60-67]. Developments on the buckling 
of unstifTened laminated composite circular cylinders may be found in references [68-77], 
In 1947, Van der Neut [78] studied the effects of eccentric stiffeners on the buckling of 
circular cylindrical shells. The work presented in [78] showed a factor of three in the 
difference between the theoretical buckling loads for internally and externally stiffened 
shells under axial compression. Baruch and Singer [79] presented work on the general 
instability of a simply supported cylindrical shell under hydrostatic pressure that was 
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analyzed by considering the distributed stiffness' of the frames and stringers separately, 
taking into account their eccentricity. Additional theoretical work on the buckling of 
isotropic cylindrical shells with eccentric stiffeners may be found in the papers by 
Hedgepeth and Hall [80], Singer et al. [81,82], Block et al. [83], and McElman et al. [84], 
Some of the first experimental work on the buckling of eccentrically stiffened cylinders 
was conducted by Card [85] and this work was compared to theoretical results by Card 
and Jones [86], Many other papers on the theoretical and experimental buckling of ec- 
centrically stiffened cylindrical shells are available in the open literature. The calculation 
of accurate buckling loads for stiffened composite cylinders is a formidable task because 
of material anisotropy, various loading and boundary conditions, skin-stiffener inter- 
action, differing moduli in tension and compression, and nonlinear behavior. 

Analysis of stiffened laminated cylindrical shells was first employed using the smeared 
stiffener approach. This type of analysis treats the eccentrically stiffened composite shell 
as an equivalent laminated cylindrical shell. A variational procedure is usually employed 
in order to obtain the results. The smeared approach was used by Simitses [87-89] and 
Jones [90,91] for the stability analysis of ring and stringer (axially) stiffened composite 
cylindrical shells. Simitses [87-89] considered the stiffened circular cylindrical shell as 
being orthotropic and reduced the strain-displacement relations to the Donnell type 
equations. Various loading conditions such as axial compression, lateral pressure, 
hydrostatic pressure, and torsion are considered for shells with clamped boundary con- 
ditions in references [87-89]. Jones' work [90,91] was presented for a circular cylindrical 
shell with multiple orthotropic layers and eccentric stiffeners under axial compression, 
lateral pressure, or a combination thereof. Classical stability theory which implies a 
membrane prebuckled state was used for the simply supported edge boundary condi- 
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tions. More recently, Reddy [37] has developed a smeared approach for axial and ring 
stiffened composite shells using the layerwise theory'. 

A new technology known as continuous filament grid stiffening has enabled the manu- 
facture of complex stiffened cylindrical shells. This cost effective process reduces the 
number of parts and fasteners since the stiffeners are integrally wound as part of the 
shell. In this study, emphasis will be upon geodesically stiffened shells produced by the 
aforementioned manufacturing process. To date, published work on the subject of 
geodesically stiffened shells is sparse. Buckling analysis of orthotropic cylindrical shells 
with eccentric spiral-type stiffeners using the smeared technique was conducted by 
Soong [92] for simply supported shells subjected to one of the following loadings: axial 
compression, hydrostatic pressure, torsion, and bending. Soong concluded that based 
on equal stiffener weight or equal strength, the spirally stiffened cylinders are about 
equal to the ring and stringer cylinders for axial compression and pure bending loads, 
but are superior in resisting torsion hydrostatic pressure loads. Meyer [93] studied an 
isotropic geodesically stiffened shell have 45 0 integrally milled out one sided stiffeners. 
This type of stiffener arrangement was used to exclude the buckling modes between hoop 
reinforcements. Meyer used a smeared approach for simply supported shells and con- 
cluded that no increases in axial critical loads were obtained for addition of internal 
pressure. Studies of isogrid composite cylindrical shells were conducted by Rehfield et 
al. [94] as well as Reddy et al. [95] extended the work to orthogrid stiffened composite 
shells. In both papers [94,95] a Donnell type theory was used for general instability, skin 
buckling, and stiffener buckling. Shaw and Simitses [96] used a smeared procedure in the 
nonlinear analysis of axially loaded laminated cylindrical shells with various in place 
transverse supports. The work in [96] includes the effects of geometric imperfections and 
lamina stacking sequence. Further work on geodesically type stiffened cylindrical shells 
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using a smeared approach may be found in references [97-99]. The smeared stiffener 
technique is effective if the cross sections of each stiffener is the same and the stiffener 
spacing is small. If the number of stiffeners is small or the spacing is large, the smeared 
stiffener analysis does not yield accurate results and usually overpredicts the buckling 
load. 

A procedure other than the smeared technique must be used for the buckling analysis, 
vibration and/or stress analysis of sparsely stiffened shells. It is desirable to treat the 
stiffeners and skin as separate structural components to determine the most accurate 
buckling or vibration mode and the local peak stresses and strains. The discrete analysis 
procedure is the only alternative to a finite element analysis to study localized effects. 
Several authors [100-103] have studied the vibration analysis of discretely stiffened cy- 
lindrical shells. Because of the relatively simple geometry of ring stiffened cylindrical 
shells, treatment of the circumferential rings as discrete elements have been considered 
in several papers [104-107]. Wang et al. [108,109] first developed a discrete analysis for 
isotropic cylindrical shells with stiffeners and then later extended the same concepts to 
composite cylindrical shells with stiffeners [110]. In the discrete analysis of [110] separate 
equations are developed for the axial stiffeners, ring stifTeners, and skin. The equations 
are coupled through interacting normal and shear loads via the application of an Airy 
stress function to the compatibility relations. Pochtman and Tugai [111] used a discrete 
analysis to study the stability of composite cylindrical shells stiffened with cross ribs. The 
development was based on the principle of minimum potential energy where the strain 
energy of the skin and the stiffeners were treated as separate quantities. Chao et al. [ 1 1 2] 
also employed the principle of minimum potential energy in the analysis of stiffened 
orthotropic foam sandwich cylindrical shells. The authors in [112] included the effects 
of transverse shear deformation in their development. Birman [113] applied a discrete 
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analysis to the divergence instability of reinforced composite cylindrical shells. The de- 
velopment consisted of solving the equations of motion in terms of displacements. The 
Dirac delta function was applied to discretely include the stiffeners' extensional, bending, 
and torsional terms in the equations of motion. Additional references on buckling of 
discretely stiffened cylindrical shells may be found in (114-117]. 

Another method of constraining stiffeners to the skin is by the application of the 
Lagrange multiplier method. Several authors used the Lagrange multiplier method in 
plate stability problems in order to satisfy boundary conditions (118-121). Al-Shareedah 
and Seireg (122] correctly predicted the transverse deflection of a pressure loaded rec- 
tangular isotropic plate with an oblique stiffener. Lagrange multipliers were used to en- 
force transverse displacement continuity between the plate and stiffener at a finite 
number of points. Phillips and Gttrdal (123] applied the same technique to the stability 
of orthotropic plates with multiple orthotropic oblique stiffeners. The Lagrange multi- 
plier method should be viable for stiffened composite circular cylindrical shells. Johnson 
and Rastogi (124] applied the Lagrange multiplier method to orthogonally stiffened 
composite cylindrical shells in order to determine the interacting loads between the 
stiffeners and the shell wall when the shells are subjected to internal pressure. No studies 
are presented in the open literature on the buckling of stiffened layerwise plates or shells 
having discrete stiffeners using an analytical method. The Lagrange multiplier method 
could easily be used to attach the stiffeners to the skin of a layerwise plate or shell. 

Finite element analysis of stiffened structures has been divided into several categories. 
The simplest yet least accurate method is to use a coarse model with lumped stiffeners. 

In the lumped stiffener method each stiffener is lumped into the plate or shell on the 
nearest element boundary. The stiffeners are assumed to be connected along the nodes 
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of the plate or shell elements as bar elements. This model introduces inconsistencies. 
The lumped method is theoretically inaccurate, as the lumped stiffener indicates a cou- 
pling along the nodes to which it is connected whereas a stiffener placed within a plate 
or shell element indicates coupling of all the nodes of the element. Further, diagonal 
stiffening is difficult to achieve with this method. A second approach is to use 
orthotropic simulation (smeared technique) of stiffened structures. This method and its 
deficiencies was discussed previously for buckling analysis of stiffened cylindrical shells. 
Another approach is the development of a special bending element where the stiffener 
stiffnesses are incorporated into the bending element at the elemental level see references 
(125-130]. This method may work well for bending, but the effects of in-plane loadings 
are not documented. Also, obtaining the skin/stiffener interaction mechanisms is diffi- 
cult to extract using this approach. The final method of modeling stiffened structures 
is by representing the stiffeners as beam, plate, or shell elements. This method provides 
the greatest accuracy, the most realistic model of skin/stiffener interaction, and conse- 
quently will be the method used in this work. 

When employing the discretely stiffened finite element approach often curved beams are 
used as reinforcing members for shells. The beam elements must have a compatible dis- 
placement pattern with that of the shell. Kohnke et al. (131] analyzed an eccentrically 
stiffened cylindrical shell by using a beam finite element with displacements compatible 
with the cylindrical shell element. Venkatesh and Rao [132] developed a laminated 
anisotropic curved beam finite element to be used in conjunction with anisotropic shell 
elements [133-135]. Bhimaraddi et al. [136] used shear deformable laminated curved beam 
elements to study stiffened laminated shells. Ferguson and Clark [137] developed a var- 
iable thickness curved beam and shell stiffening element with transverse shear deforma- 
tion for isotropic elements. Reddy and Liao [138,139] utilized degenerated 3-D beam 
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elements as stiffening members in their nonlinear analysis of composite shells. An alter- 
native approach to stiffening shell structures with beams is to approximate the stiffeners 
with the same element type used for the shell (140J. Using this procedure results in the 
introduction of a substantial number of additional nodes and nodal displacements. 

Work on postbuckling analysis of stiffened shells is sparse. Knight and Starnes et al. 
(1,2] have done some work on the postbuckling analysis of stiffened and unstiffened 
composite panels using a finite element analysis. Sandhu et al. [141] performed a finite 
element analysis of the torsional buckling and postbuckling of composite geodetic cyl- 
inders. This work concluded that joint flexibility is an important factor in the overall 
shell behavior. Hansen and Tennyson 1142] presented an overview of the development 
of a computer model for analyzing the crash response of stiffened composite fuselage 
structures. A finite element formulation was presented that supposedly can treat lami- 
nated shell buckling, large deflections, nonlinear response, and element failure. However, 
no results were presented for this work. 

The displacements, stresses, and failure analysis of shells is receiving more attention than 
in the past. Leissa and Qatu [143] applied the Ritz method to study the stresses and 
deflections in composite cantilevered shallow shells. Boitnott, Johnson, and Starnes 
[144] calculated the linear and nonlinear interlaminar stresses for pressurized composite 
cylindrical panels. The work in reference [144] also included a nonlinear failure analysis 
of pressurized composite panels. Failure was found to occur near the corners of the 
panels along the boundary of the panel. Research work on the stress distribution near 
the stiffener intersections is lacking. The layerwise theory could easily be adapted to the 
analysis of stresses near the stiffener intersections. Of particular interest may be the 
interlaminar stress at the stiffener intersections. Furthermore, the layerwise theory is a 
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quasi 3-D theory which overcomes the finite element aspect ratio problem of traditional 
3-D elements. 


1.2.3 Failure Mechanisms 

Failure analysis of stiffened composite structures is a highly complex and sparsely re- 
searched area. The failure scenarios for stiffened composite structures include: general 
instability (global structural buckling), stiffener buckling (crippling), skin buckling, and 
material failure. If the structure is designed to have postbuckling strength, then failure 
will most likely be based upon ultimate rather than buckling strength. Spier [145] con- 
ducted a failure/column buckling analysis of graphite epoxy stiffened panels using a 
mechanics of materials approach. A comparison of skin buckling, stiffener crippling, and 
structural buckling was made. Reddy et al. [95] performed an analysis based on me- 
chanics of materials in their study of isogrid and orthogrid stiffened composite circular 
shells. Their analysis considered general instability, rib (stiffener) crippling, and skin 
buckling. 

In order to determine the failure load of a stiffened structure some type of failure crite- 
rion must be applied. Two approaches to failure may be used. The mechanistic (micro- 
mechanics) failure approach deals with the failure of a composite material at the 
constituent material (fiber, matrix) level. The micromechanics approach is difficult and 
often the results are intractable except for simplistic models. The phenomenological 
(macromechanics) failure prediction is developed by treating the composite as a homo- 
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geneous material where the effects of the constituent materials are detected only as av- 
eraged composite properties. 

The mode of failure of laminated composites may be by fiber yielding, matrix yielding, 
fiber failure, delamination, or fracture. The first three failure modes depend on a con- 
stituent s strength properties. Delamination generally occurs in the form of cracks in the 
plane of the composite, resulting from manufacturing defects, low strength of resin rich 
regions, and high local stresses due to improper stacking sequence. Fracture is the result 
of preexisting voids or cracks in the constituent materials. Macroscopic failure criteria 
are based upon the tensile, compressive, and shear strengths of an individual lamina. 

A myriad of literature exists concerning failure of composite materials. A survey of 
macroscopic failure criteria applied to composite materials is presented by Sandhu [146], 
Tsai [147], Tsai and Hahn [148], and Nahas [149]. Some of the more popular failure cri- 
teria include the maximum stress criterion, maximum strain criterion, and quadratic 
polynomial criteria such those proposed by Hill [150], Azzi-Tsai [151], Chamis [152], 
Hoffman [153], and Tsai-Wu [154]. The maximum stress criterion and maximum strain 
criterion are called independent mode failure criteria and thus there is no interaction 
between modes of failure. The quadratic polynomial failure theories are mathematical in 
nature and are basically empirical curve-fitting techniques. There exists considerable 
failure mode interaction with the polynomial failure theories. The Tsai-Wu criterion is 
a tensor failure theory which is invariant under rotation of coordinates and transforms 
via known tensor transformation laws. None of the aforementioned failure criteria can 
predict the mode of failure. Hashin [155] proposed a failure criterion which considers 
four distinct failure modes - tensile and compressive fiber and matrix modes. 
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Several authors have presented some relatively simple micromechanics failure ap- 
proaches. Craddock and Zak [156] developed a theoretical model which accounts for 
large transverse stresses in the plies (laminae) and permits gradual plastic yielding of the 
matrix to failure. Sanders et al. [157] applied simple micromechanics failure models such 
as microbuckling, 'kink-band' failure, layer shear, and various interactive modes for ap- 
plication to composite aircraft design. 

Initiation of failure is often determined via the first-ply failure analysis. Cope and Pipes 
[158] conducted finite element analyses of composite spar-wingskin joints and ultimate 
strength was predicted through application of Tsai-Wu, maximum stress, and maximum 
shear failure criterion. Reddy and Pandey [159] conducted first-ply failure analyses of 
composite laminates. The maximum stress, maximum strain, Hill, Tsai-Wu, and 
HofTman failure criterion were used in their analyses. Kim and Soni [160,161] developed 
an analytical technique to predict the onset of delamination in laminated composites. 
Their work was extended by Brewer and Lagace [162] to develop a quadratic delami- 
nation criterion. This criterion is an average stress criterion which compares the calcu- 
lated out-of-plane interlaminar stresses to their related strength parameters. The 
criterion showed excellent correlation with experimental delamination initiation stresses. 
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1.3 Present Work 


The literature review presented in the previous sections indicates that analysis of stiff- 
ened composite shells is an area of extreme interest. The major emphasis of this research 
is to develop numerical techniques to study the buckling. linear, nonlinear, and failure 
behavior of geodesically stiffened circular cylindrical shells. The layerwise laminate the- 
ory of Reddy (LWTR) will be extended to stiffened circular cylindrical composite shells. 
Developments using the LWTR for shells will be applied using both a Ritz variational 
technique and a finite element approach. Application of appropriate failure criterion 
will be applied to the model in order to determine the appropriate failure scenario. 

The present study was undertaken with the following objectives: 

1. Develop a layerwise Ritz variational method with discrete stiffeners using the 
Lagrange multiplier constraint approach. Use this method to study the buckling of 
axially, ring, and geodesically stiffened cylindrical composite shells. 

2. Develop and verify a layerwise finite element algorithm for accurate prediction of 
displacements and stresses in composite plates and shells. The stiffeners are to be 

modeled as layerwise beam elements. Linear and nonlinear strain displacement re- 
lations are to be considered. 

3. Calculate the displacements, in-plane stresses, and interlaminar stresses in stiffened 
cylindrical shells with emphasis on geodesically stiffened shells when the shells are 
subjected to various loading conditions. 
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4. Apply failure criteria to study the first-ply failure of geodesically stifiened cylindrical 
composite shells. 

The governing equations of stifiened laminated shells using a layerwise theory is pre- 
sented in Chapter 2. Chapter 3 deals with the development of the Ritz variational 
method and the Lagrange multiplier constraint method. The finite element model, ele- 
ment types, numerical approach, and finite element verification problems are presented 
in Chapter 4. The results for several problems are described in Chapter 5. Conclusions 
and recommendations are presented in Chapter 6. 
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Chapter 2 

Governing Equations 


2.1 Introduction 

The development of refined shell theories for laminated composite shells has been moti- 
vated by the shortcomings of the classical lamination theory. The classical lamination 
theory (CLT) as applied to shells is based upon the KirchhofT-Love hypothesis in which 
the shear deformations are neglected. Consequently, first-order and higher order theories 
were developed to account for transverse stresses. These theories provide improved 
global response for deflections, natural frequencies, and buckling loads. However, these 
theories which are based upon a continuous and smooth displacement field do not yield 
good estimates of interlaminar stresses. Improved theories must be applied to model the 
local behavior near stiffener intersections of stiffened shells because laminate failure 
modes may depend upon the interlaminar stresses. The layerwise laminate theory of 
Reddy (LWTR), which has been shown to work well for plates, will be extended to cir- 
cular cylindrical shells. The basic equations for circular cylindrical shells using the 
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LWTR will be presented in the next section. Also, included in the development will be 
governing equations for discrete stiffeners. 


2.2 Displacements and Strains for Laminated Shells 

The LWTR is a displacement based theory in which the three-dimensional elasticity 
theory is reduced to a quasi three-dimensional laminate theory by assuming an ap- 
proximation of the displacements through the thickness. The displacement approxi- 
mation is accomplished via a layerwise approximation through each individual lamina. 
A polynomial expansion with local support (finite element approximation) is used in this 
development. Consider a laminated circular cylindrical shell with N orthotropic lamina 
having the coordinate system described in Figure 2. The displacements u, v, w at a ge- 
neric point (x, y, z) in the laminate are assumed to be of the form (see Reddy [37]) 

N+ 1 

;» i 
/V+ I 

v(x,y, z) = Vj{x,y)<tf(z) - ( 2 . 1 ) 

7-1 
<V+ 1 

M.*,y, 2) - Yj w M<y)<tJ{2) - w/fr 1 

7-1 


where N is the total number of layers (N + 1 interfaces including the surfaces) and 
Uj, Vj, wj are undetermined coefficients. The & are any continuous functions that satisfy 
the condition through the entire thickness 
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1^(0) = 0 for all y = 1,2 (jV+1) (2.2) 

In this development the summation convention will be used for repeated subscripts and 
superscripts. 

The approximation in Eq. (2.1) can be viewed as the global semi-discrete finite element 
approximations [163J of u, v, and w through the thickness. The ft denote the global 
interpolation functions, and u jt v>, and W] an: the global nodal values of u, v, and w at the 
interface locations through the thickness of the laminate. A finite element approximation 
based on the Lagrangian interpolation through the thickness can be obtained from Eq. 
(2.1). In this study a linear interpolation will be assumed and thus 

“i =(/<*> 

v, - V\'\ v 2 - - V ?\ ... pf (2.3) 

w, * W^\ w 2 = »*> = W j 2) W f > 

where W\ m fVf*> represent the values of U, V, W at the i-th node of the k-th lamina 
as displayed in Figure 3. 

The linear global interpolations are given by 

f*?" '>(*). (A-2,3 iV+ 1) 

z k^ z ^2 k+ , (A = 1, 2, ... , iV) 

where (i- 1, 2) is the local Lagrange interpolation function associated with the i-th 
node of the k-th layer as defined in Figure 3. 
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The theory presented in this study will be based upon the circular cylindrical shell ana- 
logue of the von Karman large deflection theory. This theory was applied by Donnell 
(5) and assumes that the lateral displacement w is large enough that in plane forces and 
displacements must be considered in the nonlinear form. The strain displacement re- 
lations of the Donnell/von Karman type [5,63] are 
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(2.5) 


where R is the radius of the circular cylindrical shell. The coordinated system used in 
this analysis is defined in Figure 2. Upon substitution of Eq. (2.1) into Eq. (2.5) yields 
the following relationships 
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for all i,j = 1 , 2 N + 1. 

2.3 Displacements and Strains for Laminated Beams 

The layerwise theory is extended to beams in a procedure similar to laminated shells. 
Consider a laminated beam comprised of N orthotropic lamina having a coordinate 
system described in Figure 4c. The displacements u and w at a generic point (rj, <f, £) in 
the beam are assumed to be of the form 

1 

"to, 0 - £ upi)4 ?{ o = Uf 

y - 1 

n+\ (2.7) 

h o - X 

y-i 

Here the u is the local displacement alonp the axis of the beam and w is the transverse 
displacement. In this research the out-of-plane stiffness and subsequently the out-of- 
plane displacement v is generated from the ratio of the out-of-plane beam bending mo- 
ment of inertia to the in-plane beam bending moment of inertia. See section 4.5 for a 
description of the out-of-plane stiffness generation. Torsional stiffness of the layerwise 
beams is not inherently present in the layerwise elements, but this could be included if 
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l-igure 3. Variables and interpolation ructions for the shell layers. 




an assumed displacement distribution through the thickness of the stifTeners was made. 
Including torsional stifTness would involve significantly more development for the beam 
elements and is not included in this study. The importance of torsional stifTness in 
layerwise stifTeners needs further study and will be left for future work. 


The strain relations for the stifTeners are developed in a procedure similar to that of the 
shell. The stifTeners are modeled as discrete structures and thus developments are made 
for individual stifTeners. A description of the stiffener coordinate systems is provided in 
Figure 4. Figures 1 and 5 contain illustrations of the geodesically stiffened and axial/ring 
stiffened shells respectively. The stifTeners are assumed to behave like beams. In addi- 
tion, the displacement field is assumed to be similar to that of shells. See references 
[132-137,164,165] for similar curved beam developments. The stiffener strain displace- 
ment relations of the Donnell/von Karman type are 


m 

e CC s 

Sr 


du 1 / dw \ w 

dr, + 2 V dy ) + R' 


dw 

dC 

du dw 
dC dr, 


( 2 . 8 ) 


u 

R' 


Here the radius of the stiffeners /?', is developed from vector calculus and analytical ge- 
ometry (see Figure 6) and is given by the following relation 


R' = 


R 


(2.9) 


where 


Stiffener pitch = 2wb 


( 2 . 10 ) 
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Here b is the pitch parameter as shown in Figure 6. For ring stiffeners b = 0 and 
R' = R and for axial stiffeners b = oo and thus R' - oo. 


Upon substitution. of Eq. (2.7) into Eq. (2.8) yields the following layerwise following re 
lations for the stiffeners 






c ii - w t ' 


dp 

d C 


( 2 . 11 ) 


d<t> 1 d\Vf , u, , 
c «= u ‘— 


for all i,j 1, 2, ..., number of beam interfaces. 


2.4 Variational Formulation for Laminated Shells 

The principle of virtual displacements will be applied to the shell and stiffeners sepa- 
rately. For the shell, the principle of virtual displacements can be used to derive a con- 
sistent set of differential equations composed of N constant thickness l.nUn, The 
Principle of Minimum Potential Energy dn may be expressed in variational form as 

<511 = S U + <5 V =* 0 (2 12) 

Here SU is the virtual strain energy (virtual work done by the internal stresses) and S V 
represents the variation of the potential of the applied forces. The minimum potential 

energy statement for the shell in terms of stresses and virtual strains caused by virtual 
displacements may be expressed as 
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Figure 4. Coordinate systems for the stilTcners: (a) axial stiffeners; 







an 



{a ] Sc ] + <j 2 Sc 2 + o 3 dc 3 + <t 4 <5c 4 + o 5 6c 5 + o 6 Sc 6 }dzdQ + 6 V = 0 (2.13) 


where a u a 2 , a it <*♦, a s , and a 6 are the shell stresses, <5e,, Sc 2 , <5e 3 , <5c 4 , da and Sc 6 are the 
virtual strains in the shell coordinates, h is the total shell thickness, Q denotes the total 
shell area at the midplane, and once again 6 V represents the variation of the potential 
of the applied forces. The variation of the total potential energy in terms of the stress 
resultants, displacements, and virtual displacements is obtained by substituting the strain 
displacement equations (2.6) and integrating through the shell thickness. The variation 
of the potential energy then takes the following form 
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where the stress resultants, ML, M'l Q\, Q\ J, 0, K\, and and the variation of the po- 
tential of the applied forces, 5 V , take the following form 


Governing Equations 


32 



(«= 1 , 2 , 6 ) 


K = o a <t>‘dz , Ml = J_ 2 oj^dz , 

T T 

h i 

(Qh (? 2 - Ql) = ( ct 5 > <Ji, 03) ~~ dz 

~ 

± 

(A^,A^)= f ((J 5 , o A )<p‘dz 

y —h 
2 

± 

SV= ~J P^wdQ. ~ J J 2 {N n du n + N s 5u s )dzdQ. 

Q ~2~ 


(2.15a) 


-£L 

f fa/ A dSw _ a 5(5w. a \ 

- J 0 L( m " ~ +M, ~ + 


for a linear prebuckling analysis <5 V reduces to 
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For the potential energy of external loads p is the applied pressure, N„ and N, are the 
applied in-plane normal and tangential forces respectively, and M „ , M„ and Q are the 
applied edge normal moment, tangential moment, and shear force respectively. For a 
buckling analysis N\, N 2 , and A'is are the axial, lateral, and shear external forces respec- 
tively acting on the shell membrane. 


The cylindrical shell is assumed to be laminated of orthotropic layers with the principal 
material coordinates of each layer oriented arbitrarily with respect to the shell axis. The 
layer constitutive equations referred to the shell coordinates are given as 
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(2.16a) 



where Q are the components of the orthotropic stiffness matrix, 
in terms of displacements is given by the following expression 
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The stress resultants 
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for all i,j,k,l - 1, 2, ..., N + 1 and where 
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Note that Dfy, Dfy, Di$, and Z)#' are symmetric in their subscripts and superscripts. 



e/c. 


(2.19) 


The coefficients with a single bar over them are not symmetric with respect to the 
superscripts. The variational statement in terms of displacements is provided in Ap- 
pendix A. 
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2.5 Variational Formulation for Laminated Beams 

The variation of the potential energy for the beam along the beam length, L„ may be 
expressed in the following manner 


<5n 



a m bz rm + a U 6c Zi + 


SV 


( 2 . 20 ) 


The three dimensional constitutive equations for an anisotropic body are reduced to that 
of a one-dimensional body by eliminating the normal stress the in-plane shear stress 
a vt , and the transverse shear stress <T tf . Similar procedures for the modeling of laminated 
composite beams were employed by Bhimaraddi and Chandrashekhara [164] and more 
recently by Kassegne (165]. The stresses <7 f{ , a^, and are eliminated, but the strains 
C{{» and cjc are not eliminated. For a laminated beam the constitutive relations re- 
duce to the following form 



( 2 . 21 ) 


where the components of the reduced orthotropic beam stiffness matrix, C$, are ex- 
pressed in terms of the original orthotropic stiffness terms, Q, and are expressed in the 
following form 
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Finally, the nonlinear variational statement for laminated beams in terms of displace- 
ments may be expressed as 
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for all i,j,k,l = 1, 2 f number of beam interfaces* 


Here Qj, represents the in-plane area of the beam elements and where 
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Here h is the beam height and integration is made through the height of the beam dz. 

Note that Z%, B'Jp, /?$, and B'jj! 1 are symmetric in their subscripts and superscripts. 
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The coefficients with a single bar over them are not symmetric with respect to the 
superscripts. 


The potential of the external forces for a beam is given as 
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ter all i,j = 1, 2, ..., number of beam interfaces. 


(2.26) 


where L, is the length of the beam and F, is the force acting along the length of the 
stiffener. 
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2.6 Failure Equations 


The various failure criteria were discussed in section 1.2.3. In this study, a macrome- 
chanics based first-ply failure analysis will be conducted for some select cases. As dis- 
cussed previously there are many macromechanics based failure criteria. The failure 
analysis involves calculating the stresses and strains at a point in the structure and then 
applying the selected criterion. These criterion include the the experimentally deter- 
mined macroscopic material strength data. In this research work, the Tsai-Wu failure 
theory is used as the working failure criterion. The Tsai-Wu criterion was selected be- 
cause of its general character. The Tsai-Wu criterion has three distinct advantages: (1) 
invariance under coordinate rotation; (2) transformations are made via known tensor 
transformations; and (3) there exists symmetry of properties similar to those of the 
stiffnesses and compliances. Therefore, the Tsai-Wu criterion was selected for this work. 

The Tsai-Wu criterion is given by the following expression 


F i°i + F^OfOj £ 1 i, j * 1 6 (2.27) 


Here a, are the stress components and F, and are the strength terms. The strength 
tensor terms may be expressed as 
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( 2 . 28 ) 


All other strength tensor components are zero. Here <7|, <72* 03 are the normal stress 
components, <t 4 , <t 5 , <t 6 are the shear stress components, X T (X C ), Y T {Y C ), Z T {Z C ) are 
the lamina normal tensile (compressive) strengths in the x, y, z directions respectively, 
and R, S, T are the shear strengths in the yz, xz, and xy planes respectively. The values 
for Xr, X c , Y t , Y c , Z t , Zc, R, S, and T will be given later in this research work. 
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Chapter 3 

Ritz Buckling Method 


3.1 Introduction 

In this chapter a method is developed to study the buckling of stiffened cylindrical 
composite shells with discrete stifFeners using a closed form analytical solution. The 
stifTeners are directly attached to the shell where the components of the displacements 
between the shell skin and the stiffeners is accomplished via the application of the 
Lagrange multiplier method. Many of the equations developed in Chapter 2 are appli- 
cable to the Ritz buckling method derived here, but some simplifications are also incor- 
porated. In lieu of layerwise beams, the Euler-Bernoulli beam theory is used in 
developing the discrete stiffeners. The method developed in this chapter is applicable to 
cross-ply and some quasi-isotropic shell layups. In this study, simply supported edge 
boundary conditions will be assumed in order to apply the closed form solutions. 
Equations (2.1)-(2.6) and (2.12)-(2.19) are applicable for the layerwise shell used in this 
approach. The Euler-Bernoulli beam theory is presented in the next section. 
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3.2 Euler-Bernoulli Beam Stiffeners 


The governing equations for the Euler-Bernoulli beams are developed in a procedure 
similar to that for layerwise beams of sections 2.3 and 2.5. The stiffeners are modeled 
as discrete structures and thus the development of a generic stiffening element may be 
applied to axial, ring, or geodesic stiffeners. The displacements for a stiffener using the 
Euler-Bernoulli beam theory are given by the following relation 
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'=1,2 T 


(3.1) 


where u' s and w } are the displacements for each of the T stiffeners. Here u is the local 
displacement of the stiffeners. For ring stiffeners the u's are replaced by v's. A de- 
scription of the stiffener coordinate system is provided in Figure 4. 

The Euler-Bernoulli strains for the stiffeners are developed from the displacements and 
are given as 
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The definition of R ' was developed in section 2.3. Here £ is the distance from the 
stiffener centroid to the reference surface. The uniaxial stress-strain equation for the 
stiffeners is 
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where <r' is the stress in an individual stiffener, £' is the modulus of the stiffener, and t 
is the strain in an individual stiffener. 


The variation of the potential energy for the stiffeners may be expressed as 

«5n = SU‘ s + 6V* S /= 1, 2 T 


(3.4) 


The variation of the strain energy SU‘ s for an individual stiffener may be expressed as 



(3.5) 


where 


• /i moment of inertia of an individual stiffener about the reference surface (ie., ( 
= 0) and thus /' = I' xe + (Ci) 2 ^j. 


• Cl - distance from the stiffener centroid to the reference surface. 
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• Lj = length of an individual stiffener. 

The variation of the external forces for the Euler-Bernoulli beam stiffeners is developed 
from the potential energy statement and is expressed as 


, .j [ L ’\ -a dw‘ s d5w‘ s ] 

iv ‘~ l \ + 


(3.6) 


where F, is the force acting along the length of the stiffener. 


3.3 Lagrange Multiplier Method 


The procedure used in applying the fundamental mathematical principles of Lagrange 
multipliers is described briefly in this section. The Lagrange multiplier method will be 
used to constrain the discrete stiffeners to the shell surface. Proof of the validity of the 
Lagrange multiplier method is provided in references [118, 119]. Al-Shareedah and 
Seireg [120-122] successfully applied the Lagrange multiplier method to stiffened plates. 
Phillips and GUrdal [123] applied this method in the stability analysis of orthotropic 
plates with multiple orthotropic geodesic stiffeners under in-plane loadings. A de- 
scription of the Lagrange multiplier method follows. 

Let it be required to minimize a function of I variables ... , xj 

A*\,* x,) (3.7) 
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where the x s are not independent but are bound together by J independent constraint 
functions hj such that 


^1 C^l • *2’ "■ » X j ) = 0 

h 2 (x l ,x 2 ... ,x,)^0 

(3.8) 


hj{x u x 2 ... , Jf/) = 0 


Lagrange's method of simultaneously minimizing f and satisfying the constraint func- 
tions is to minimize a composite function L called the Lagrangian defined as 
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where x denotes the vector of variables x\ f xi , ... , jr/anda is the vector of undetermined 
Lagrange multipliers «i,a 2 , ... ,a y The necessary conditions for/fx) to be a minimum 
while simultaneously satisfying the constraints a are 
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Differentiation of Eq. (3.9) results in the following set of equations 
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Equation (3.11) is a set of (I -+- J) simultaneous equations for (I + J) unknowns 
xand a The Lagrange multipliers may be used for variational type problems where the 
\ariational statement is considered a functional similar to_/(.v). The Lagrange multipliers 
may be used to constrain the displacements and or rotations of a discretely stiffened 
structure. 


3.4 Stiffened Shell System 


In this study, the stiffeners are treated as discrete structures and are attached to the skin 
using the Lagrange multiplier constraint technique discussed in the previous section. 
See Figure 7 for a representation of a discretely stiffened shell showing the constraint 
points. The variation of the total potential energy for the stiffened cylindrical shell may 
be expressed as 


<3n — <5(/ SHELL + ^STIFF + ^^SHELL + ^^STIFF (3-12) 

The virtual work for the shell was derived and is shown in Eqs. (2.12)-(2.19) and in Ap- 
pendix A. The Euler-Bemoulli beam strain energy and potential energy was shown in 
Eqs. (3.5) and (3.6) respectively. For the buckling analysis used here only the linear 
portion of the virtual work statement is applicable and thus is given as 
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(3.13) 
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for i,j 1, 2, .... N + I and for k - 1, 2, ..., T (number of stiffeners). 

The definitions of M\, A/j, Mi, Q[, 0, 0, K\, ft, If h N 2 , and N 6 are provided in 
Eqs. (2.15) and (2.17). The development presented here is for specially orthotropic cyl- 
inders where C l6 = C 26 - C* * C 4S - 0. However, when the values of 

Cl6 ’ ^ 26 ’ ***• and C4J are sma11 such as the case for certain quasi-isotropic materials 
then the layerwise Ritz method should provide reasonable results. This is demonstrated 
in the results presented in Chapter 5. 


For a buckling analysis we have 


N x - iN x 
N 2 - xn 2 

V 6 =■ an 6 

where X is the minimum buckling eigenvalue. 


(3.14) 


We have 


* " 1» ^2 — 0, N& * 0, buckling under axial compressive load. 
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• A'i = 0, N 2 = — 1, Ai = 0, buckling under lateral pressure. 


3.5 Buckling Solutions and Equations 

The Ritz method will be employed for this buckling analysis. The results may be com- 
pared to other solution procedures like those described by Jones [90,91], Reddy [37], and 
linear finite element methods. The global buckling solution of stiffened circular cylin- 
drical shells consists of solving Eq. (3.12) such that 

<5n — 0 (3.15) 

In order to solve the equations using a Ritz variational approach, a solution must be 
assumed. In this study, simply supported edge conditions will be studied. The edge 
boundary conditions for the shell skin and the individual stiffeners are (see Figure 2) 

iv = v = 0 at jc = 0, L (3.16) 

The following solutions of the Navier form which satisfy the boundary conditions are 
assumed: 

Shell Skin 

u, - Ur cos(flt^r) cos tPj) 1 - 1 , 2 , ... ,{N + 1) 

v t = V?* sin(a wr x) sin^) ; a m L, /?„ = -£- (3.17) 

"i- sin(a m jr) cos(/? n y) 

where (//"", Vp", and Wp" amplitudes are to be determined for each mode (m, n). 
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Axial and Geodesic Stiffeners 


u l s = U‘ k cos(a^) ; * k = -SS- , i-l,2 T 

w'= W lk sin(ot*jr) 


( 3 . 18 ) 


Ring Stiffeners 

«i- ^sin(^) ; A/-—-, f-1,2. F 

w'= tffcos(/V) 


(3.19) 


where f/J* ( F'O and IFf* ( amplitudes are to be determined for each stiffener mode k 
(1). For axial stiffeners the r\ axis is along the x direction and for ring stiffeners the rj axis 
is along the y direction. 


Substituting Eqs. (3. 1 7)-(3. 1 9) into the buckling expression, Eq. (3.13), and integrating 
over the regions ( L c = 2nR) below 


<50 = f (<5 ^shell + & F SHELL ) ^ + f (^^stiff + ^stiff) 

J a J o 

= f f (<5t/sHELL+ ^SHELI.) ^ + f (^^STIFF + ^ ^STIFf) dlj 

J 0 J 0 “0 


(3.20) 


yields the following expression 
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+ (*>sW + *>U : + p- ^ + x ^ + -£- BJ + (3 ' 21 

+ ^(^a m 2 + A^,, 2 ) Wj mn }6 W™ 

+ ZT{^(**‘ l < , -«*T- f ‘ 1 ‘'<)< 

+ [^W< - 4<£(y --%■*?< ) 

For ring stiffeners replace (/]* by K(', ^'* by If"', and a* by ( - /?,). 

3.6 Constraint Equations 

In this part of the study the displacements u and w of the shell and stiffeners will be 
constrained (v and w for ring stiffeners) to attach the stiffeners to the shell skin. The 
compatibility of displacements is imposed at a finite number of points along each 
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stiffener (see Figure 7). The compatibility of displacement equations may be expressed 
as 

Axial and Geodesic Stiffener Constraint Equations 

e j P - u, - ui = £//"" cos(a m x Jp ) cos (p„y Jp ) - U> k cos(a k ij jp ) = 0 
fjp - w i ~ - Wf" s in(a^) cos(P„y jp ) - W Jk sin(a k t} Jp ) = 0 

i-l,2. (W+l) ( 3 - 22 ) 

,r 

P = 1,2, 

Ring Stiffener Constraint Equations 

- KT" sinCa^) sin(/?^) - if sin(/?^) = 0 

hj P = = W™ sin(a m jr ;> ) cosC/?^) - wf cos(^,tf Jp ) = 0 

/-1,2. ((V + 1) ( 3 - 23 ) 

;-i,2,...,r 

P=l,2 /> 

Here the subscript p represents the number of constraint points for P total constraint 
points along the stiffeners. A set of Lagrange multipliers may be developed where one 
multiplier is required for each constraint condition. The equations take the form 

1 - y Jp e Jp + nj/ jp + <f> jpgjp + + Jp h Jp (3.24) 

where y p , n p , <f> Pf and \j/ p are the Lagrange multipliers and L is the set of Lagrange mul- 
tiplier equations. The system of Lagrange multipliers L may be added to the potential 
energy FI to form a variational functional 1 such that 

/ = n + L (3.25) 
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In order to minimize the functional, the variation of I may be set equal to 0 so that 

SI - <5n + SL = 0 (3.26) 

The minimum potential energy <5n has been defined previously. The variation of the set 
of Lagrange multipliers is 

SL = 5yj p e jp + y Jp 6U™ cosfa,^) coitfitfjJ - y Jp SUi k cos(oc k rt Jp ) 

+ SHj j/jp + n Jp S W™ sin(a m jTyp) cos^^) - n Jp 6 w{ k sin (a kVjp ) 

+ &4> Jp gjp + 4> jp S V™ sin ( am x Jp ) sin (P„y jp ) - <t> Jp 5 V* sin {P nVjp ) (3 ' 2?) 

+ 5 *lP h Jp + M W ™ sin (««*y,) cosC/?^) - ij/ Jp S W J ! sin(P n rjj p ) 

From Eq. (3.26) <5/» 0 so the individual components must therefore be 0 and hence 


6U™ - 0 

Sy JP - o 

6V? n = 0 
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II 

O 

6W™ - 0 

&4>jp = o 
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II 

V 

c* 

II 

o 

SW{ k - 0 


5V{‘ = o 


On 

II 

o 



Eq. (3.28) yields 11 sets of equations which may be solved. If no ring stiffeners are 
present then 9 sets of equations must be solved. Similarly, if no axial or geodesic 
stiffeners are present then 9 sets of equations must be solved. Carrying out the oper- 
ations of Eq. (3.28) yields the following sets of equations. 
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(3.29a) 


i u -r 5 - { K. + d'Lp? + d‘L) ur - xMrt + d‘L) yp" 

~ “”( ~R D '' 1 + ®>J “ ®?j) H 7'"| + >7p cos Kr*y p ) cos(/!^, f i = 
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+ -fir- (x k Up + j + A^hj ! H^*} - >iy, sinfa,^,) - 0 

I 1 •’f + A -Jjf - + Stf “ ^JP sin ^fljr> “ 0 


0 


(3.296) 


(3.29c) 


(3.290 


(3.29c) 


(3.29 j/) 
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(3.29 g) 


+ 


R' (Pi V “ + R > ) + - i, jp cos(fi jp) = 0 


6y Jp : U r cos(a m jfy p ) cos(/?^ p ) - if cos (a k rt Jp ) - 0 

(3.29/i) 

Sfi Jp : K n sinfot^) co$(P„y Jp ) - W J 5 k sin(a^) = 0 

(3.290 

Hjp. K mn sin(a m x Jp ) sin(/?^ p ) - ff sin(fa Jp ) - 0 

(3.29 i) 

S +Jr W ‘ mn sin («m*jp) cos(/?^) - w{ 1 cosy fl Jp ) - 0 

0.29k) 


3.7 Shell/Stiffener Load Distribution 


3.7.1 Introduction 

The appropriate prebuckling load distribtion as applied to the shell (skin) and the 
stiffeners is essential for proper analysis of the stiffened shell structure. In the smeared 
buckling analyses the skin and stiffener properties are averaged to form and equivalent 
structure and therefore no prebuckling load distribution is necessary. However, for dis- 
crete structures the proper distribution of loads applied to the skin and the individual 
stiffeners must be obtained. The approach taken here for the prebuckling load distrib- 
ution involves using classical lamination theory (CLT) and a smeared load distribution 
whereby the stiffeners are considered to be smeared for the sake of calculating the ap- 
propriate prebuckling load distributions only similar to the approach used by Phillips 
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and Gurdal [123] in their study of stiffened plates. The stiffeners will be assumed to carry 
only forces along the length of the stiffener. The development of the constitutive 
equations presented here follows the approach of Jones [166], 


3.7.2 Shell Constitutive Relations 


Using CLT for the shell/skin, the stress strain relations in the principal material coordi- 
nates for a 2D laminate are 



Qn 

Q\2 

0 


Q\2 o 
0.22 0 
0 ( 266 . 



(3.30) 


Here the Q tj are the reduced stiffnesses. If the coordinate system is transformed then the 
stress-strain relations are given by 



2,. 

(2.2 

Q\ 6 

Q\2 

Q 22 

Q 26 

Q 16 

Q 26 

Q 66 



(3.31) 


The Qtj are the transformed reduced stiffness terms. Using the Kirchhoff-Love hypoth- 
esis the displacements u and v within a laminate may be expressed as 


u 

v 


* dw° 

U — 2 ■ 


0 

V — 2 


dx 

dy 


(3.32) 
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Here u°, v°, w° are the displacements at the midsurface of the shell, 
strains may be expressed as 


Vxy = 


£ 

— in 

du° d* w° 


dx 

** Z <3x 2 

= 

dv __ 

dy 

dv° w° d 2 w° 

% K 2 dy 2 

du 


du° . dv° „ 

dy 

+ & 

" ay + Bx - 21 ' 


5 2 o 
W 

dxdy 


or 




+ z 



where the middle surface strains are 


O'] 


r du° > 

dx 

d 

- { 

dv ° - . a*® 

h 

f - \ 

1 9y * R + Sx I 

Y°xy 


% 
V 


and the middle surface curvatures 



-- 


r 

d 2 u> g 

a* 2 

a2 o 
0 W 

dy 2 

2 _aV_ 

dxdy 

v. s 


The subsequent 


(3.33) 


(3.34) 


(3.35) 


(3.36) 
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Substituting the strain relations (3.34) into the stress-strain relations (3.31) the stresses 
in the k' h layer of the shell may be expressed as 




2 12 

2,6 

2,2 

Qn 

2 26 

2l6 

2 26 

2 66 



(3.37) 


The resultant forces and moments acting on the shell laminate may be obtained by in- 
tegration of the stresses in each layer (lamina) through the laminate thickness. The re- 
sulting expressions are 



where 2 * and r*_ 1 are defined in Figure 8. Upon integration through the laminate the 
following expressions result 



BU 

B s n 

Bn 

B s 22 

B \ 6 

B s 2 6 
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(3.40) 



where 


(4 4 °«) - (3.41) 

w —ft 

2 

The superscript s denotes that these are shell(skin) constitutive relations. 


3.7.3 Axial Stiffener Constitutive Relations 


The effective stiffener spacing l x over which the influence a particular axial stiffener has 
upon the structure is described in Figure 9 and is given by the expression 


2ttR 


(3.42) 


where N a is the total number of axial stiffeners. The linear strain displacement re- 
lationship for axial stiffeners is 


du aV 

™ a <2 * 

dx dx 2 


(3.43) 


where t% is the smeared axial stiffener strain and tf and are the respective displace- 
ments of the axial stiffeners. The stress strain relationship is given as 
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= Ec - E 

^a L x t-'a 


du 1 
dx 


— z ■ 


-v2 a 

<7 W 


<3* 


( 3 . 44 ) 


where <r? is the smeared axial stiffener stress and E a is the modulus of the axial stiffeners. 
The total force in the stiffener is calculated by integrating the stress distribution in the 
stiffener over the area, A a , and is expressed as 



= E a A a c a x + E a A£ a K x 


( 3 . 45 ) 


The force resultant, N x , is then calculated by distributing the total stiffener force F; over 
the stiffener spacing l x and is written as 




( 3 . 46 ) 


The moment resultant Ms, is calculated by multiplying the force resultant Ns, by the 
distance z a from the neutral axis to the centroid of the stiffener at which the force acts. 

Here it is assumed that the forces act at the centroid of the stiffeners. The moment re- 
sultant is then 



l x 


j £ x 4" » K x 

‘x l x 


( 3 . 47 ) 


where the moment of inertia l a is 
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h = 4 + ^a^a 


( 3 . 48 ) 


The set of constitutive relations for the axial stiffeners becomes 


N x = A \\ £ °x + B^K x 

K - + 


where 


Z) 


a 

11 


A 

Bn 


E a A a 

lx 


E a A a^a 

lx 



EJg 

lx 


(3.49) 


(3.50) 


3.7.4 Ring Stiffener Constitutive Relations 


The effective stiffener spacing l y over which the influence a particular ring stiffener has 
upon the structure is described in Figure 10 and is defined by the relation 



(3.51) 


where L is the length of the cylinder and N c is the number of ring stiffeners. The linear 
strain displacement relationship for ring stiffeners is 
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where c c y is the smeared ring stiffener strain and V s and w c are the respective displacements 
of the ring stiffeners. The stress strain relationship is the given as 


C 



-j c c c 

CV W _ Q W 

Sy R z Sy ! 


( 3 . 53 ) 


where a c y is the smeared ring stiffener stress and E c is the ring stiffener modulus. The 
total force in the stiffener is calculated by integrating the stress distribution in the 
stiffener over the area, A c , and is expressed as 



EgAfZy T EgAgZgKy 


( 3 . 54 ) 


The force resultant, N y , is then calculated by distributing the total stiffener force F y over 
the stiffener spacing l y and is written as 



( 3 . 55 ) 


The moment resultant M c y% is calculated by multiplying the force resultant N c y , by the 
distance z e from the neutral axis to the centroid of the stiffener at which the force acts. 
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Here it is assumed that the forces act at the centroid of the stiffeners, 
sultant is then 




o Eclc 

I c y + / K y 

y l y 


where the moment of inertia I c is 



The set of constitutive relations for the axial stiffeners becomes 


Aj - A 22 e° + B c 22 K y 
My ~ + D 2 2 Ky 


where 


*22 


B c r 


'22 


D c 


22 




b 


The moment re- 


(3.56) 


(3.57) 


(3.58) 


(3.59) 
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3.7.5 Geodesic Stiffener Constitutive Relations 


Described in Figure 1 1 is the coordinate system for the geodesically stiffened shell and 
the definition of a single geodesic cell. The effective circumferential length /, is the cir- 
cumference divided by the number of cells per shell circumference such that 


where N a is the number of geodesic cells. The cell length is given by L t . The strains in 
stiffener coordinates are 



(3.61) 


where s, is strain along the stiffener axis, e t is strain transverse to the stiffener axis, 
is the shear strain, and u* and w* are the geodesic stiffener displacements. See Figure 4 
for the stiffener coordinate system. 

The stifieners 1 and 2 of Figure 1 1 are oriented at angles 0 and - <f> respectively from 
the y axis of the shell. Since the stiffeners are not aligned with the shell coordinate sys- 
tem, the values in the local coordinate system aligned with the stiffener axes must be 
transformed to the shell coordinates. The strains in the local coordinate system aligned 
with the stiffener can be determined using the strain transformation relations by Jones 
[166]. These are 
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(3.62) 



s 1 4> 

C 2 (f) 

S<t>C(f) 

C 2 (j) 

s 2 <t> 

— 50C0 

—2 s4>c4> 

2s<t)c<t) 

C 2 0 — S 2 4> 



where 


s<t> = sin <t> 

c<{> = cos 4> (3.63) 


The axial strain component for stiffener #1 of the unit cell at angle + <f> yields 




* 2 2 

sin 4>c x + cos <t>c y + sin <f> cos <j>y ^ 
2 2 

= s 4>c x + c <t>t y + s^c^y^ 


(3.64) 


The axial strain component for stiffener #2 of the unit cell at angle - </> yields 


h = sin 2 ( — 4>)c x + cos 2 ( — <fr)z y + sin( — <f>) cos( — 0)y^, 
= s 2 <j>c x + c 2 <f>c y — s^ctfrYxy 


(3.65) 


The stiffener axial forces F\ and F 2 may be found by multiplying the strains e, and c 2 by 
the stiffener elastic modulus E t and integrating over the area of an individual stiffener 
area A t . This results in the following expressions 


F i = (s 2 <t>c x + c 2 <j>c y + sfftc^YxyjdAg 

A t 

F 2 = £ S J (s 2 <f>c x + c 2 <t>c y - s^c^YxyjdAg 

A. 


(3.66) 
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Solving for the force F x along the circumferential edge involves resolving F { and F 2 into 
components along the x direction as indicated in Figure II. 


F\ = /•*! sin 4> 4- F 2 sin 4> 


= E„ 


J ( s 3 <j>c x + c 2 <t>stf>c y + s 2 <t>c<t>y xy )dA g 


+ J (s 3 <j>c x + c 2 <t>s4>Cy — s 2 (f>c<f>y 

A t 

~ 2 ^rJ (s 3 <t>e x + c 2 <t>s<j>c y )dA g 




(3.67) 


The force resultant, Nf, distributed over the circumferential length is calculated by dis- 
tributing the force F x over the length / r Distributing the force ftover/, and carrying 
out the integration of Eq. (3.67) yields (constant cross sectional area) 


N* 


2E g A g 3< 2 E x A t , 

s 4>z x -\ C 2 <t>Stj>Z 


/. 


8 

2E„A e z 


L 


8^8 ,3, . 2 Wi 2. . 

S <f>K x H ; C 4>s4>K v 


(3.68) 


4 


1. 


where for a constant cross sectional area 



A similar procedure is used to calculate the force resultant Nf along the lengthwise edge 
and produces the following result 
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N g 
1 y 


2 E x A g , 

— 7 s'<f)c 4 >c x + 


-g 


^g^g 3 , 

J C (pC y 

^8 


+ 


2 EgdgZg 2 
— 7-^ r <t>c<t>K x + 


■"g 


^ E g^g Z g 3 , 

— 7 c <P K v 


(3.70) 




To determine the shear force resultant, the forces tangential to either the circumferential 
or lengthwise edges will produce eigenvalue results. The tangential force components 
along the edge are F, cos 0 in the positive y direction and F 2 cos 0 in the negative y di- 
rection. Substituting for F\ and F 2 of Eq. (3.66) yields 


F* 

r xy 


= £, cos 4> — F 2 cos 0 


= E 


J (•S 2 0c0e jt + c 3 0e y + s<f>c 2 <l>y xy )dA i 

A i 

- J (s 2 0c0e x + c 3 <f>c y - s<t>c 2 4>y xy )dA l 

A t 

" 2E g\ { c2 <t>s4>Vxy)<1A g 


(3.71) 


Distributing the force F „ over the length l t and integrating Eq. (3.71) yields the shear 
force resultant 


N xy ” — } c 0s0y^, + 7 c 2 0s0i c^, ( 3 . 72 ) 

'g ‘g 

The loads acting upon the stiffeners are ofTset by an amount z r Thus, moment result- 
ants are introduced into the problem. The resultant moment Mi may be calculated by 
multiplying Eq. (3.67) by z and distributing over l t which yields 
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M g x = 


zf* 2E g r 

l ~ J fax + c <l>s<i>c y ]zdA g 


2E C 

tg 


2E ? A g z g 2E g A g z g 

s <j>t x H ; c 050 C V 


2EJ, 


4 


s 


'S‘g 3 


4 


i <t>K x + 


2E b I. 


g‘g 2 




C (f)S<j)K v 


where 


4 “ J 2 ^g * 4f + 


Similar procedures will produce the moment resultants A// and A/1,, These 
sultants are 


itr* _ 2E t A fy 2, ■ 2£ rV* 3 

/ «* <pc<pc x + - c <f>e v 

u s L g r 

2E g I. 2 22L/_ , 

+ ~ 7 S 2 fa<t>K x + ~ - C 3 <f>K y 

L g L g y 


Mt 2 W 2 , . 2 Va 2j 

™ i c <f>s<f>y xy -\ C 3 <!>s4>ty 


The constitutive equations for geodesic stiffeners are given as 



A*\ ^f 2 0 

A*, 0 

0 


22 

0 Al 


66 



*?. B* n 

*f, *? 2 

0 0 


0 

0 

*i 6 



(3.73) 


(3.74) 
moment re- 


(3.75) 


(3.76) 
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where 



Bfi 0 

Bl B f 2 0 


>i\ 

0 


& 



ofi Of: 

Z)f, £>22 

0 0 



(3.77) 


A g 

^ll 


A g 


12 


A g 


21 


^22 

Al 


fif, = 


2E g A g 3 

j 0 

l g 

2E g Ag 2 
— - C 2 4>S<f> 

l g 

2E g A g 2 
5</>C0 

L g 

2E,A e 


^g 

2E„A. 


1. 


g 


12 


B lf ; 

^6 


4 

2 V* 2* A 

C 050 


/, 


B%2 

b * 66 


2E gAgZg 2 , , 

T 5 <pC<p 

L g 

2E gAgZg 3 
T C 0 

2 , , 

' C 050 


/„ 


■C 3 0 

ZJf, « 

3 1 

. 5 <p 

l g 

• C 2 050 

Df 2 = 

2 E I 

LC *g k g 2 t g 

C 050 
( g 

^-s 3 0 

$ 

II 

2E g I g 2 

l g 

— C 2 050 

D22 = 

2EgIg 3. 
, c <p 


(3.78) 


'g 


3.7.6 Skin/Stiffener System Constitutive Equations 

Using the principle of superposition the force and moment resultants for the stiffened 
structure are 
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(3.79) 


'Vx = A* + N* x + N* 

N y = Ny + N c y + X 8 y 

A xv = Ax, + A* 

The stiffness terms (A, B, D) are then 

^11 = + 4 a u + 4 * t 

^12 = + 4* 2 

4\6 = ^*6 
^26 = ^26 

^22 = ^22 + ^22 + 4*2 
466 = ^66 + 4*$ 

&U = + 5*, + 5f, 

5,2 = 5f 2 + 5f 2 
5, 6 = 5[ 6 


Af x = M x + M x + M* 
M y = m; + m; + m* 

A/xy = + A/f, 


-#26 = ^ 2 6 

5 22 = 5 2 2 + ^22 + #22 
566 = 5^ 6 + 5 f 6 

A, = £u + Z>n + 0fi 

D,2 = D s n + Z)f 2 (3.80) 

5*16 = D\ 6 
5*26 = ^26 

5*22 = 5*22 + 5*22 + 5>22 
5*66 = 5>g 6 + D* 6 


It is assumed that the skin and stiffeners have identical strains. The constitutive 
equations for the shell (skin)/ stiffeners are: 




3.7.7 Loading Conditions 


Four different loading conditions will be considered for the solution of shell/ stiffener 
prebuckling load distributions; axial compression, pressure loading, applied shear 
(torsion), and applied end shortenings. In this study only buckling due to axial com- 
pression and/or pressure loading will be considered. The other two loading conditions, 
torsion and applied end shortening, are included for completeness and can be included 
in future work including calculating the prebuckling load distributions for a finite ele- 
ment analysis. Load distributions for combinations of the above loading conditions such 
as for example axial compression and pressure loading may be obtained via the super- 
position of two or more load conditions. 


3.7.7. 1 Case 1 - Axial Compression (Applied Nx) 

For this case it will be assumed that the prebuckling N y =N xy = 0. Further, due to the 
offset loads acting at the stiffeners' centroid, a resultant moment is developed. An ap- 
plied N x may be reacted by either axial or geodesic stiffeners. Thus, the resulting mo- 
ment is 


M x ™ 2 a Nl + z.N* 


(3.82) 


The resulting offset moment in terms of strains is 

M x = ZaiJnZ + B^k x ) + I g (A* u c x + Af 2 c° + Bf x K x + Bf 2 tc y ) (3.83) 

Substituting Eq. (3.83) into Eq. (3.81) reduces the constitutive equations to 
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The set of equations, (3.84), will be solved using a Gauss elimination procedure for 
CS, Cy, Y°xy, *x, Kxy» The solution for the skin force resultants is then 

= AnC x + A^Cy + A^y^ + l K x + B\2 K y + &\6 K y 
Ny = ^]2 c x + ^22 c y + ^26 Y°xy + B\2 K x + &22 K y + #26 K y (3.85) 

^Xy ~ d l6 E X 4 - A^ly + 4 - B^ 6 K X + B^ 6 Ky + + BfaKxy 

The resultant force in the axial stiffener is 


K - A* u i° x +B a u K x (3.86) 

and the force in each axial stiffener is 

K - l x K (3.87) 

The resultant force in the geodesic stiffeners is 

N* - A* n t° x + A* n t° y + Bf ]Kx + Bf 2 K y (3.88) 

and the force on the the geodesic stiffeners is then 
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( 3 . 89 ) 



2 s<p 


It is assumed that che ring stiffeners do not carry any loads when the shell is subjected 
to axial compression so F c y = 0. 


3.7. 7.2 Case 2 - Pressure Loading (Applied Ny) 

For pressure analysis the value of AT, - -pR where p is the external pressure. If internal 
pressure is applied then N y = pR. For this case it will be assumed that the prebuckling 
Nx=N *y = °- Due ^ the offset loads applied to the stifTeners this results in an applied 
moment resultant. This moment resultant is created because of the ofTset force result- 
ants acting on the ring or geodesic stiffeners acting at the centroid and may be expressed 


M y = Z C N^ + Z g Ny 


( 3 . 90 ) 


The resulting offset moment in terms of strains is 

My - zlA C n t 0 y + B C n Ky) + Z^A^ + A^ty + 5 ^ + ^Ky) 

Substituting Eq. (3.91) into Eq. (3.81) reduces the constitutive equations to 
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(3.92) 



The set of equations, (3.92), will be solved using a Gauss elimination procedure for 
c 5. £ y. Vxy, k x , Ky, K xy . The solution for the skin force resultants may be found from Eq. 
(3.85). The resultant force in the ring stiffener is 

Ny = A 22 c y + B 22 K y ( 3 . 93 ) 


and the force in the ring stiffeners acting at the ring stiffener centroid 


~ ^yi^22 c y + &22 K y) 


(3.94) 


The resultant force in the geodesic stiffeners is 


‘Yk — ^2l c x + ^22 £ y + + &22 K y 


(3.95) 


and the force on the the geodesic stiffeners is then 


F\ 


- ^2 - 



2 c<f> 


(3.96) 


The axial stiffeners carry no circumferential load and hence F$ * 0. 
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3.7.7.3 Case 3 - Shear Load (Applied Nxy) 


In this case a state of applied shear loading, N xy , exists. For this case it will be assumed 

that the prebuckling N x = N y = 0. The applied offset moment resisted by the geodesic 
stiffeners is then 


M xy - VV (3.97) 

The resulting offset moment in terms of strains is 

M xy = + BiiKxy) (3.98) 

Substituting Eq. (3.98) into Eq. (3.81) reduces the constitutive equations to 



The solutions of Eq. (3.99) for c$, c®, y£,, k x , tc y , and k v may be substituted into Eq. 
(3.85) to find the skin resultant forces. The solution for the skin force resultants may 
be found from Eq. (3.85). The load distribution in the geodesic stiffeners may be found 
from the following expressions 

^xy = = l g{ A l(>y 0 xy + ^>(> K xy) (3.100) 
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(3.101) 


Fy - -F, 


pg 
1 x y 

2 aj) 


The axial and ring stifTeners carry no shear loads and hence F; = 0 and F c y - 0. 


3.7.7 .4 Case 4 - Applied End Shortening 

In this case the strain, c?, must be calculated from the applied end displacements A such 
that 



(3.102) 


The solution will involve solving for N x as an unknown rather than cj. It is assumed that 
in the prebuckling state N y , N v , M y , and are all zero. The offset moment will be 
given by Eq. (3.83). Substituting Eq. (3.83) into Eq. (3.81) and rearranging to solve for 
N x as an unknown with ej known yields the following equation 
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(3.103) 



The set of equations. (3.103) will be yield ft * y} r , K „ K „ K „. The solutions and the 

known Cl may be substituted into Eqs. (3.85X3.89) to calculate the appropriate forces 
in the skin and stiffeners. 


3.8 Governing Equations and Final Form 

The set equations of 3.29 represent 11 sets of equations which may be written in matrix 
form as 

CXii*i|](A) - (3.104) 

The nonzero elements of the stiffness and mass matrices of Eq. (3.104) are provided in 
Appendix B. The stiffness matrix of Eq. (3.104) is sparse and thus in this study the 
eigenvalue problem will be solved in terms of the shell (skin) displacements U. V, and 
w. The Lagrange multipliers. y„ and are eliminated from the problem by 
expressing these terms as functions oTU, V, and W. After reducing the 11 x 11 set of 
equations the following eigenvalue problem results 

C 5 n3 iS n l CS 13 ] fupl '[0] [0] [0] I 

[S| 2 ] [S v ] CS 23 3 <{V} \ - x [0] [0] [0] ){V) } (3.105) 

[^ 3 ] [5 23 ] [S 33 ]_ [0] [0] [Af 33 ] 
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where 


[s,,] = [^,] + w, 8 ][^ 8 r , [w,4r l [^,] 

[Sij] = CA^ia] 

[5 13 ] = c^3] + c^. 8 ][A'4. 8 r , [A;5][A 9>5 r 1 cA 9> 3] 

[5 21 ] = [^21 ] 

C 5 223 = [* 22 ] + [^2. 1o]C^5.1o3 -1 [^<563[^1 0 ,6] -1 [^10.2] 

[•^ 23 ] ~ C^233 + [^2, lo3[^«, lo3 C^673[^ll,73 * C-^l 1. a3 

[5 31 3 = [^.3 + CA 3> 93[^5 1 93" i [^J[A 8 ,J‘ 1 [A Si1 ] 

M = [^32] + [^3J«]C^.t,3" , [A' 76 3[^0,6]" , [Ar i 0.23 

[S 333 - ^333 + [^3.93C^5.93" , [^S53[^.s3" 1 CA 9(3 ] 

+ [^.,«. 1 ,]" , C^7 7 3C^„.73- , [^„.33 

[A/ 33 ] - [^333 + CM 3>9 ][A/ 5i9 ]- , [iVf 55 ][A/ 9 5 ]- , [A/, i 3 ] 

+ [A/3,i.]CA/ 7 ,„r , CM 77 ][A/ 11 , 7 r , [Af n( 3] (3,I06) 

The solutions of Eq. (3.105) yields the eigenvalues X„„ for each mode M, N and the 
minimum eigenvalue is the critical buckling load. 

In this development the out-of-plane and torsional stiffnesses of the beams were neg- 
lected. However, both the transverse and the in-plane displacements were constrained 
between the shell skin and the stiffeners. In many similar analyses only the transverse 
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displacements were constrained [118-123] and reasonable results were achieved. More- 
over, this development is just an initial study of the buckling of stiffened shells using the 
layerwise theory and an analytical approach. The results shown in Chapter 5 reveal that 
this method does work well. Hence, including the out-of-plane and torsional stiffness 
of the stiffeners may in fact overstiffen the structure as developed in this chapter. This 
could yield poor results when compared with other analyses. The majority of this re- 
search involves the development and the use of the layerwise finite element method de- 
scribed in Chapter 4. Consequently, it was decided to concentrate more upon the finite 
element method and analyses rather then to include every minute detail into the analyt- 
ical buckling approach. Including the out-of-plane and torsional stiffnesses of the 
stiffeners can be included in the future. 
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Chapter 4 

Finite Element Formulation 


4.1 Introduction 

Displacement based finite element models developed from the governing equations pre- 
sented in Chapter 2 are derived for both layerwise shell and beam elements. The 
layerwise theory reduces the equations of three-dimensional elasticity to a quasi three- 
dimensional laminate theory by assuming a layerwise approximation of the displace- 
ments through the thickness. Consequently, the strains are different in different layers. 
The interlaminar stresses ( a xt , a yn o„) will be calculated using information from the in- 
plane stresses calculated from the finite element solution and by using an approximate 
technique to integrate the equilibrium equations. The variational statements of Appen- 
dix A and Eq. (2.23) are used in the development of the finite element models. In addi- 
tion, the derivation of the direct stiffness and tangent stiffness matrices are presented. 
The finite element method for plates and shells is discussed in Refs. (32, 41-59, 125-142, 
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163, 167-170). Several example problems are included in this chapter to verify the finite 
element program. 


4.2 Layerwise Shell Finite Element Formulation 

The generalized displacements (u Jt v Jt vv,) for the shell elements are expressed over each 
element as a linear combination of the two-dimensional interpolation functions i p n and 
the nodal values (u ; ", vf, wf) as follows 


/vos 

( U J> v p w j) = v j> w j)'i ,n 


n m 1 


(4.1) 


where NDS is the total number of nodes per element. Substituting Eq. (4.1) into the 
variational statement of Appendix A yields the shell finite element model. A geometric 
description of the finite element shell model is shown in Figure 12a. The elemental finite 
model for layerwise shells may be expressed as 

A I1 A 12 

a 21 a 22 
j Lrfnn rsmn 

A 31 A 32 

where K^\ n , ATpy , etc. are the element submatrices provided in Appendix C, (U), (V), and 
{W} are the column vectors of nodal displacements, {F„}, {F,}, and {F,} are the column 
vectors containing the boundary and force contributions. 
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Figure 12. Geometry of the finite element model: a) shell element; b) beam element. 
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Three types or finite elements are developed in the finite element program. These are 

linear (4 nodes), Serendipity (8 nodes), and Lagrange (9 nodes) elements. Figure 13 

shows these three types of finite elements. In this study the isoparametric rectangular 

master elements ace used. The interpolation functions for these elements are listed be- 
low. 


Linear Element ( NDS = 4 ) 

" T (1 “ W ~ *) 

n ) = -“(1 + 0(1 - r ,) 

*) - jd + Od + 17) 

*> " T (I " + v) 

Serendipity Element [ NDS =» r ) 

At. "> - yd - fld - ,x - 1 - {- ,) 

>1) - yd + «1 -,)(-!+{-,) 
At. 1) - yd + i)d + u)( - 1 + { + ,) 
At. v) - |d - {xi + ix - 1 - { + ,) 

^ 5 ({. 9) - yd - { ! XI - 7 ) 

At. 1) - y(l + {)(! - A 
At. 1 ) - yd - { ! xi + n) 

At. 9) - y(l - {)(! - A 
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Lagrange Element ( NDS = 9 ) 


As, >;) = y(S 2 - S)(>f 2 - »r) 

As, 17) = y(S 2 + S)(if 2 - v) 

As, «r) = y(S 2 + S)(if 2 + 1?) 

As, if) = y(S 2 - «v a + v) 

As, *) - yd - S 2 )(if 2 - If) ( 4 . 5 ) 

As, if) - y(S 2 + 0(1 - i/ 2 ) 

As, >f) - y(l - S 2 )(if 2 + if) 

AS, if) = y(S 2 - S)(l - if 2 ) 

As. 17) = (i - Ad - if 2 ) 


4.3 Layerwise Beam Finite Element Formulation 


The generalized displacements (uj, wj) for the beam elements are expressed over each el- 
ement as a linear combination of the one-dimensional interpolation functions i// n and the 
nodal values (uf, wf) as follows 


NDS 

(«y, Wj) = £<* W jW 

n m 1 


( 4 . 6 ) 
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2 


(- 1 .- 1 ) 
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(- 1 , 1 ) 

( 0 , 1 ) 

( 1 . 1 ) 
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(- 1 , 0 ) 
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9 
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$ 

1 

( 0 , 0 ) 
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Figure 13. 


Nod, numbering art f„ Li.enr, Smart, art Ugrnnge .bell finile elemen... 


Finite Element Formulation 


90 



where NDS is the total number of nodes per element. Substituting Eq. (4.6) into the 
variational statement of Eq. (2.23) yields the shell finite element model. The elemental 
finite model for layerwise shells may be expressed as 

hmn l^mn 
A tl A 1 2 

Lfnn Lmn 
_ A 21 a 22 

s s 

where ATI", etc. are the element submatrices provided in Appendix C, {U} and {W} 
are the column vectors of nodal displacements, {F„} and {F c } are the column vectors 
containing the boundary and force contributions. See Figure 12b for a description of the 
layerwise beam element geometry. 

Three types of finite elements are developed in the finite element program. These are 
linear (2 nodes), quadratic (3 nodes), and cubic (4 nodes) elements. Figure 14 shows 
these three types of finite elements. In this study the Lagrange family of master elements 
are used. The interpolation functions for these elements are listed below. 

Linear Element ( NDS — 2 ) 


w 

{W} 


sw 


(4.7) 


= yd - 0 

^ = -*-(1 + i ) 


(4.8) 
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Quadratic Element ( NDS = 3 ) 


+' - - yid - £) 
t 2 - (1 + «(l - *) 

<A 3 = -±-«l + i) 

Cubic Element ( NDS = 4 ) 

- -jy(‘ - «y + «Xy - <) 

- -fyO + «Xi - {X-f - {) 

* J - -fg-ti + {Xi - {)(y + « 

- - -jy(y + {Xy - {)(! + {) 


(4.9) 


(4.10) 


4.4 Assembly and Nonlinear Analysis 


The layerwise shell and beam elements are assembled directly into the global stiffness 
matrix which yields the following sets of equations 

[K(A)]{A} - {F) (4.11) 

such that 

WW)] - (CAJ + [^(A)]) (4.12) 

where [AT(A)] is the assembled global stiffness matrix, [tfj and [*m.(A)] are the linear 
and geometrically nonlinear parts of the global stiffness matrix respectively, {A} is the 
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CUBIC ELEMENT 



QUADRATIC ELEMENT 


Figure 14. Node numbering and coordinates for Linear, Quadratic, and Cubic beam finite elements. 
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column vector of nodal displacements, and { F] is the column vector of boundary or ap- 
plied force conditions. The system of equations represented by Eq. (4.11) can be solved 
directly for linear problems. However, for geometrically nonlinear problems the stiffness 
matrix is a nonlinear function of the unknown solution and must be solved iteratively. 

The method selected here for solving geometrically nonlinear problems of stiffened plates 
and shells is the Newton-Raphson method. In the Newton-Raphson iteration method, 
the basic equations for the residual vector {/?} is given by the expression 

{/?} s [tf(A)]{A} - {F) - 0 (4.13) 


Assuming that the solution is known at the r ,h iteration, the residual vector {/?} is ex- 
panded in a Taylor series about {A r }, 


0 * {*}, + 


5{A} 


{<5A}, + 


- (WA’)]{Al - {F» + 


(4.14) 


where the tangent stiffness matrix [AT 7 } is given by 

M-[w] 


The tangent stiffness matrices for layerwise shell and beam elements are presented in 
Appendix C. The final equation to be solved for the increment of the solution {<5A r } is 

= {/?} (4. 16) 


The total displacement at the (r + 1) M iteration is given by 
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K +, l = (at + {^ 


( 4 . 17 ) 


The convergence criteria used to determine when the iterative solution stops is 




5S EPS 


(4.18) 


If applied displacements rather than applied forces are specified for a particular problem 
then the total loads at the (r 4* 1) ,A iteration is given by 

fO - {n+{«n (4i9) 

The subsequent convergence criteria is then 



where I is the number of nodes at which the applied displacements are not specified. 

A geometric explanation of the Newton- Raphson technique for a one-dimensional 
problem is provided in Figure 15. 
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4.5 Beam Element Stiffness Transformations 

Beam elements may be oriented at any arbitrary angle a from the x axis of the shell as 
shown in Figure 16. The translational degrees of freedom of the beam u' and vv' are re- 
lated to the displacements in the shell coordinates u, v, and w by the following vector 
relationships 


u' = wcosa -I- vsina 

v' = —u sin a 4- vcosa (4.21) 

w' = w 


The stiffness transformation relation may be expressed as (see Cook [170] ) 


ra = m’iwjm 


(4.22) 


where [AJ is the transformed beam stiffness matrix, [7] is the transformation matrix, 
and [A' J is the stiffness matrix to be transformed. For example, the beam node local to 
global coordinate transformation matrix is 


C r, 


( 3 * 3 ) 


] 


ca sa 0 
- sa cot 0 
0 0 1 


(4.23) 


where 


sa = sin a 
ca = cos a 


(4.24) 


The stiffness matrix to be transformed is 
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(3x3)3 ~ 


0 


0 

A 

K\ 


12 


«. 0 

^21 0 K' 


22 


(4.25) 


where F? u represents the in-plane stiffness, fC l2 and fC 2l are the shear stiffnesses, fC 22 
is the transverse stiffness, and A^n is the out-of-plane stiffness. 

In this research the stiffeners are thin (0.2") and an approximation of the out-of-plane 
stiffness is made based upon the ratios of the out-of-plane moment of inertia, I„, with 
the in-plane moment of inertia, I<. See Figure 12b for a description of the beam element 
and the geometry. The out-of-plane stiffness is then developed as 


where / { = -yy- and I n = 


hbl 

12 ‘ 


Carrying out the matrix multiplications as defined by Eq. (4.22) yields the following 
transformed stiffness matrix 


Z*, 


(3*3)- 


(A^,,c 2 a + A^,,s 2 a) (A",,ca.ia - £,,ca.sa) A7 12 ca 
(/f n ca.ra - £ n c<ua) (^ n j 2 a + £ n c 2 ac) 

^2i c “ fC 2] sa /C 22 


(4.27) 


Each node in the beam will transfer in the same manner as that presented in Eqs. 
(4.21-4.27). 
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u = u cos a + v sin a 


v 7 = -u sin a 
w' = w 


Figure 1 6 . Representation of the beam displacements (if, v\ W) 
mcnts (u, v, w). 


v cos a 


shell transformation shell displace- 
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4.6 Interlaminar Stress Calculation 


In this study an approximate technique is used to integrate the equilibrium equations 
by using the in-plane stress information provided by the finite element solution. The 
technique presented by Chaudhuri and Seide [171] is extended here to quadrilateral 
isoparametric elements. The work of Reference [171] is derived for interlaminar shear 
stresses and was also adopted in the work by Barbero and Reddy [172]. The authors of 
[172] obtained the transverse shear stresses using derivatives of in-plane stresses that 
were calculated by differentiating the interpolation functions of a finite element approx- 
imation based on a generalized laminated plate theory. The work presented in [171,172] 
will be extended to both transverse shear stresses (o x „ a yt ) and the transverse normal 
stresses a tt for layerwise shell elements. Additional references on calculating 
interlaminar stresses may be found in [173-175]. 

In this study the interlaminar stress distribution through each layer is approximated with 
a quadratic function requiring 3N equations for each of the interlaminar stresses 
(&xit Qyi\ Gzz ) where N is the number of layers; N equations are used to satisfy the N 
average shear stresses on each layer. Two equations are used either to impose vanishing 
shear stresses at the top and bottom surfaces of the shell or for the interlaminar normal 
stress o u = p 0 on the surfaces. If there is no applied p 9 on one or both surfaces then 
°ti - 0 . Then, (N - 1) equations are employed to satisfy continuity of the stresses as 
the interfaces between layers. The remaining (N - 1) equations are used to compute the 

dflxi dOy t do,, , . 

jump in — , — , or — — at each interface. 

The average stresses on each layer are computed from the constitutive equations and the 
displacement field obtained in the finite element analysis. 
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The equilibrium equations for a cylindrical shell are 


da xx 

dx 


+ 


1 


da 


xy 


da 


yy 


R By 

da 


R By 
d t 

1 7 (cT ~ 


+ 


xy 


d _°xi 

dz 

do 


yz 


dx 


dw 

dx 


1 


do 


yz 


)- 


da 


= 0 

= 0 

a yy 1 d 

R R dy 


dz 


R dy 


XI 1 d_ 

dx R dy 


dw 
R dy 


) 


, dw , 

( ^aT> 


w 


d , a xy d 
dx R dy 


) + ■ 


da 


22 


dz 


+ Po = Q 


( 4 . 28 ) 


The equilibrium equations for a flat plate reduce to 


da 


XX 


da 


xy 


da. 


dx 

• + 

dy 

+ 

* 

l 

dOyy 

t 

daxy 

i dtJy * i 

dy 

T 

dx 

dz 

d °X2 

+ 

dOyj 

da ZI 

+ -t £L + j 


dx 


dy 


dz 


+ A> = 0 


( 4 . 29 ) 


Here the equilibrium equations are used to compute f” , ■ ■ - — anc j di rec tly 

dz dz dz 

from the finite element approximation. The components of the stresses and their deriv- 
atives are computed from the constitutive equations for each layer. The procedure re- 
quires computation of the second derivatives of the displacements {u jt v Jt w,) are presented 
in Appendix D. 
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4.7 Finite Element Verification Analyses 


4.7.1 Introduction 


Several representative problems are analyzed using the previously derived layerwise ele- 
ments and solution procedures. Some of the problems have analytical solutions or the 
analysis has been presented in the open literature by using different finite element mod- 
els. Comparisons of the present results with published solutions, where available, pro- 
vide a check for the accuracy and applicability of the layerwise elements developed for 
this research work. Although additional analyses were performed only a selected group 
of representative sample problems are presented here. 


4.7.2 Unstiffened Plates and Shells 

1. Orth otropic Clamped Cylindrical Shell 

A comparison of the center deflection of an orthotropic clamped cylindrical shell sub- 
jected to internal pressure as shown in Figure 17 is presented in Table 1. A comparison 
is made with the finite element solutions presented by Reddy (167] and Rao [176] and the 
analytical solution presented by Timoshenko and Woinowsky-Krieger [61]. A 2x2 mesh 
of 9 node layerwise shell elements was used for this analysis. The layerwise finite ele- 
ment results yield a good correlation with the published solutions. 

2. Cylindri cal Shell Roof Under the Action of Self- Weight 
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R - 20 inches 
L - 20 inches 
p 0 - 6.4 /* psi 


Figure 17. A damped cylindrical (hell subjected to internal pressure. 
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0.25 



This problem as shown in Figure 18 was solved using conventional techniques by 
Scordelis and Lo [177]. This particular cylindrical shell problem has been used fre- 
quently to assess shell finite element performance [45, 178-181). The shell is supported 
on rigid diaphrams and is loaded by its own weight. This is a test case of the application 
of the full process to a shell in which bending action is severe due to the supports re- 
straining deflection at the ends. In Reference [178], the authors showed that using fully 
reduced integration yields more rapid convergence and better accuracy than selectively 
reduced integration on the transverse shear terms only. The results presented in Figures 
19 and 20 are for layerwise shell elements with fully reduced integration and using lxl, 
2x2, and 3x3 meshes 9 node elements. The layerwise elements produce excellent corre- 
lation with the analytical solution of Reference [177]. 

3. Center Deflection of a Simply Supported Orthotropic Cylindrical Roof 

The geometry, boundary conditions, material properties, and results for the simply sup- 
ported orthotropic cylindrical roof is shown in Figure 21. The nonlinear results were 
developed from the Newton- Raphson procedure discussed in Section 4.5. The results 
are in good agreement with the results presented by Palmerio [182]. 

4. Center Deflection of a Simply Supported 0/90 Cylindrical Roof 

The geometry, boundary conditions, material properties, and results for this problem are 
found in Figure 22. The Newton- Raphson technique was employed to acquire the 
nonlinear results. The results agree well with the work by Palmerio [182]. 
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Figure 18 . An isotropic cylindrical shell roof under self-weight. 
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Figure 19. Axial deflection at the support of an isotropic cylindrical shell roof under self-weight. 
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Figure 20. Transverse deflection at the support of an isotroj 







P/1000 (lb) 


P 


X 



\ 


Boundary Conditions : 

u - 0 at x - 0 
v « 0 at y ■ 0 
v ■ w « 0 at x - L. -L 
u - w - 0 at y - S, -S 



Figure 22. Simply supported |0/90| cylindrical roof. 
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5. Stress Analysis or a Simply Supported fO 90/01 Plate 


The plate being studied is loaded with a uniformly distributed transverse load. The plate 
boundary conditions and material properties are described in Figure 23. The plate is 
simply supported on all four sides and due to symmetry only a quarter of the plate is 
modeled. A 4x4 mesh of 9 node elements are used in this analysis. The through the 
thickness distribution of the inplane normal stress, a xx , for an aspect ratio of a/h = 10, 
is shown in Figure 24. The stresses were computed at the Gauss point x = y = 0.0528a. 
Figures 25 and 26 contain similar plots of the interlaminar shear stresses a yt and a xl , 
respectively. In Figure 25, a yt is computed at the point x = 0.0528a and y = 0.9472a. 
In Figure 26, a xz is computed at the point x = 0.9472a and y = 0.0528a. In these plots, 
dashed lines represent stresses obtained from the constitutive equations, while the 
smooth solid line represents the stress distribution obtained using the equilibrium 
equations as developed in Section 4.6. Stresses obtained using the LWTR, the FSDT 
(first order shear deformation theory) are also compared in these plots. The transverse 
normal stress, a cl obtained from both the constitutive and equilibrium equations is 
shown in Figure 27. The transverse normal stress is obtained at the Gauss point x = y 
= 0.0528a. Modeling each layer in the composite plate as several layers may serve to 
increase the agreement between the LWTR equilibrium and the LWTR constitutive re- 
sults. 

4.7.3 Beam Structures 

1. Cantilever Beam Subjected to an End Load 

The beam dimensions and properties used in this analysis are found in Figure 28. The 
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[0/90/0] Simply Supported Plate 

E,» 1 8.5 x TO 6 psi 
E 2 = E 3 * 1 .64 x 1 0 6 psi 
6 12 * ^13 * 0*87 x 1 0 6 psi 
G 23 = 0.54 x 1 0 6 psi 

v 12« v 13 “0.30 
v 23 * 0.49 


Figure 23. Simply supported |0/90/0| square plate subjected to a uniformly distributed load. 
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Figure 24. Through-tbe-thickness distribution of the in-plane normal stress for a simply supported, 
|0/90/0| laminated square plate under uniform load, (a/h “ 10). 
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Figure 26. Through-the-thickncss distribution of the transverse shear stress o n for a simply supported, 
|0/V0/0| laminated square plate under uniform load, (a/h ■ 10). 
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Figure 27. 
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linear results for this analysis are presented in Table 2 and good correlation exists be- 
tween the finite element and the classical solutions. 

2. Cantilever Beam Subjected to a Uniform Load 

The beam dimensions and properties used for this example problem are shown in Figure 
28. The linear results for this analysis are presented in Table 2. The linear finite element 
results compare well with the results from classical beam theory. The large deflection 
analysis of the cantilever beam subjected to a uniform load is presented in Figure 29. 
These results compare well with the results presented by Liao [183]. 


4.7.4 Stiffened Structures 

1. Analysis of a Stiffened Plate with Eccentric Stiffeners 

This problem shown in Figure 30 was analyzed by Liao [183]. Analyses were made using 
both 2 and 4 layers for the plate and 2 beam layers. Reduced integration was used for 
the transverse shear stiffness terms. The results of this problem are shown in Table 3. 
The LWTR finite element results compare well with the finite element results presented 
by Liao and with classical beam theory. 

2. Cantilever Plate with Symmetric Stiffeners 

The geometry, material properties, and loading condition of the cantilever stiffened plate 
are shown in Figure 31. The results obtained are compared with the finite element re- 
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E - 1 .2 x 1 (J psi 
v -0.2 
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Table 2. Linear Results for a Clamped Beam Subjected to an Applied End Load and to a Uniformly 
Distributed Load (E “ 1.2 x 10*, » m 0.2, L ■ 10 in.) 


Loadng 

Conditio n 


LWTR (FEA) OasscaJ Basn Theory 

(Irtfies) (inches) 


End Load 


P - 1 «». -0.334694 


-0.333333 


Distibutad Load 


q- 1 byin. -0.125691 


-0.125000 
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Figure 29. Large deflection of a cantilever beam under a uniform load (E ■ 1 2 * 10* psi, * - 0.2, L 
■ 10 in.). 
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suits obtained by Liao [183J and are provided in Table 4. A good comparison exists 
between the two finite element solutions. 

3. A Square Plate Resting on Elastic Edge Beams and Supported at the Corners, Sub- 
jected to a Uniformly Distributed Load 

Figure 32 shows the plate geometry and material properties. The same problem was 
solved by Timoshenko [61] who assumed that the elastic edge beams are of zero torsional 
rigidity. Liao [183] also solved this problem using a finite element technique. The results 
are obtained with a 2x2 mesh of 9-node shell elements and four 3-node elements are 
displayed in Table 5. Results compare favorably with those of Liao. 

This concludes the finite element verification analyses. The next chapter deals with the 
buckling and stress analysis of stiffened composite shells, with emphasis on geodesically 
stiffened cylindrical composite shells. 
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E • 3 x 10 7 pH 

v-0.3 

Arti - 52 ^ 2 
»n»- 152.41 h 4 
P-2lbs. 


Figure 30. Cantilever stiffened plate subjected to on end load. 
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Table 3. Transverse Deflection of an Eccentrically Stiffened Plate. 
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Figure 31. Cantilever stiffened plate with symmetric stiffeners. 
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Table 4. Transverse Deflection of a Cantilever Stiffened Plate with Symmetric Stiffeners. 


** — *- . 

Mon * 


4x4 Shall 

12 Beam 


Liao (1987)* LWTR (FEA) * 

(g|0 (cm) 

-0.18103 -0.18482 (2 layers) 


* Transverse Deflection at the Loaded Point 
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Ptat* 19.53" x 19.53" x 0.2" 

Edge beams : 0.5" width, 1" depth 
E« 10 7 psi, v > 0.25 

P-1 pci 



0 


Figure 32. A square plate resting on elastic edge beams. 
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Table 5. Transverse Deflection of an ElasticaJIy Supported Plate Subjected to a U niformly Distributed 
Load* 


Mesh Liao (1987) LWTR (FEA) 

(Indies) (inches) 


Timoshenko 

(inches) 


2x2 Shell -0.095957 

8 Beam 


-0.097258 (2 layers) -0.120290 

-0.102152 (4 layers) 
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Chapter 5 
Results 


5.1 Ritz Buckling Results 

In order to validate the LWTR for discretely stiffened shells, some numerical results are 
needed. Before the LWTR is used to generate results for discretely stiffened shells, nu- 
merical results for certain known configurations are generated and compared with the 
published solutions. A comparison of buckling results for unstiffened circular cylindrical 
shells with simply supported boundary conditions are presented in Table 6. The results 
of the LWTR are compared with an analytical solution first for an isotropic aluminum 
circular cylinder subjected to axial compression. Also, in Table 6 a comparison of the 
nondimensional buckling load for the LWTR, first-order shear deformation theory 
(FSDT), third-order shear deformation theory (HSDT), and classical lamination theory 
(CLT) for various cross-ply circular cylinders is presented. Results for the FSDT, 
HSDT, and CLT are from the paper by Khdeir, Reddy, and Frederick [185]. Results for 
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the LWTR analysis compare well for isotropic and cross-ply unstiffened circular cylin- 
drical shells as can be observed from Table 1, 

Next, a study of a- geodesically stiffened quasi-isotropic [ -45/45/90/0] s plate subjected 
to axial compression was made. The geodesically stiffened plate was presented in the 
work by Phillips and Gurdal [123]. They used a Lagrange multiplier approach to 
discretely attach beam-like orthotropic (isotropic) stiffeners to the plate and classical 
lamination theory to model the plate skin. Also, they conducted a finite element analysis 
of the plate using 9-node combined membrane and quadrilateral elements for both the 
skm and stiffeners. The stiffened plate was comprised of symmetric stiffeners with 
stiffener heights of 0.5 , 0.75", 1.0", and 1.25". A description of the plate geometry and 
the applied loads is shown in Figure 33. The boundary conditions were chosen so that 

u = w = 0 at y = 0, L y 

V - w = 0 at X = 0, L x (5.1) 

( L x - 80in., L y = 28in.) 

The Ritz solutions which satisfy these boundary conditions are 

Plate 

u, = U? n cos (a^) sin (p„y) i * 1,2, ... ,(N + 1) 

v, = V? n sin^) cos^) ; a m = -22L, = (5 2) 

w i= W” n sinCa^) sin]/?^) 

where Ur, Vr, and IV r amplitudes are to be determined for each mode (m, n). 
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Geodesic Stiffeners 


u s V s cosfa^) ; a k = -^- ) i=l,2,...,T 

, , lk Ls (5.3) 

= W * sin(a^) 

The procedure used to develop the buckling equations is exactly the same as that de- 
scribed for shells. Material properties for the plate and stiffeners are presented in Table 
7. Using the same geometry and loading a comparison of the buckling results is pre- 
sented in [123] and the LWTR for discretely stiffened plates is shown in Table 8. The 
LWTR compares well with the finite element method, and seems to produce better re- 
sults than the Lagrange Multiplier Method (LMM) of Phillips and Gilrdal, especially 
at lower stiffener heights. One reason for the difference between the buckling loads ob- 
tained from the finite element method and the LWTR discrete method could be that for 
smaller stiffener heights the finite element (plate) stiffeners have more of an effect on the 

skin than does the LWTR. This results in higher finite element buckling loads at lower 
stiffener heights. 

Next, buckling analyses of quasi-isotropic [ -45/45/90/0]* circular cylindrical shells with 
eccentric axial, ring, and geodesic stiffeners were conducted. The material properties 
used for these analyses are the same as those found in Table 7. A nominal shell radius 
of 85 and a shell length of 100” were selected. A shell thickness of 0.2” and a stiffener 
thickness of 0.2” were used. These same dimensions were use by Gilrdal and Gendron 
[186] in their design optimization analysis of geodesically stiffened shells. Comparisons 
of the discrete LWTR approach with the smeared CLT method proposed by Jones 

[90,91] and the smeared LWTR method of Reddy [37] were made for axial and ring 
stiffeners. 
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Table 6. 


LnstifTened Buckling Results. 

Buckling of a Simply Supported Unstiffened Aluminum Circular Cylinder 
Subjected to Axial Compression 


Buckling Load Roark (184] 

Buckling Load LWTR 


<lbs./inch) 

(Ibs./inch) 

Error 

6408.29 

6106.10 

-4.95% 


E - lOx 10*psi 
v ■ 0.30 

thickness - 0.30 * 


Buckling of Unstiffened Composite Cylinders 

Comparison of layerwtse Lammate Theory of Reddy (LWTR) with Theories of Ref. [IBS] 

- 2 
N L 

(L/R - 1. R/h-10), N- 

100h3E 2 


Larranation 

Theory 

N Simply Supported 

0 0 
0790 

LWTR 

0.1 S23 


HSOT 

0.1687 


FSDT 

0.1670 


CLT 

0.1817 

o 

o 

o 

0/90 /0 

LWTR 

0.2814 


HSOT 

0.2794 


FSDT 

0.2813 


CLT 

0.4186 

o°/w°/ ... 

LWTR 

0.2728 

10 layers 

HSOT 

0.2896 


FSDT 

0.2898 


CLT 

0.339S 


HSOT - Higher Order Shear Deformation Theory (third-order) 
FSDT - First-order Shear Deformation Theory 
CLT - Classical Lamination Theory 

E/E 2 - 40, G 12 -G u -0.6Ej, ^j-O.SEj, v, 2 - 0.25 
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Table 7. Material Properties Used in the Stiffened Buckling and Finite Element Analyses. 


Properties 

Values 

Ei 

18.5 x 10 6 psi 

e 2 

1.64 x 10 6 psi 

E 3 

1.64 x 10 6 psi 

G\2 

0.87 x 10°psi 

G 13 

0.87 x 10 6 psi 

G 23 

0.54 x 10 6 psi 

Vl2 

0.30 

Vl3 

0.30 

V 23 

0.49 

X T 

1 82.8 x 1 0 3 psi 

Xc 

210.5 x 10 3 psi 

Yt — Zt 

27.2 x 1 0 3 p«i 

Y c = Z c 

17.6 x 10 3 psi 

R 

13.5 x 10 3 psi 

S = T 

21.75 x 10 3 psi 

ply thickness 

0.005 in./ply 

ZsTIFF 

18.5 x 10 6 psi 

stifTener thickness 

0.20 in. 


Results 





Table 8. Analysis of [ — 45/45/90/0]5 0.2" Thick Plate with Geodesic Stiffeners Subjected to Axial 
Compression N, (Lx ” 80", Ly ” 28", 12 Stiffeners). 


Stiffener 

Haight 

(inches) 

... 

Taetbad FEA 
Hmt 

(lbsjlnch) 

Obfc/hdi) 

N lwnn 

LWTSD 

(beyinch) 


n lwtto 


N FEM 

0.5 

573 

342 

512 

0.596 

0.894 

0.75 

705 

613 

656 

0.870 

0.930 

1.0 

748 

708 

713 

0.946 

0.953 

1.25 

783 

i 

743 

784 

0.949 

1.001 


Finite Element Budding Load 
CLT Lagranga hMtipBer Budding Load 
LWTR Discrete Bucking Load 
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Results for quasi-isotropic [ —45/45/90/0]^ shells with external (eccentric) axial stiffeners 
subjected to axial compression are shown in Tables 9 and 10 and in Figure 34. The re- 
sults in Tables 9 and 10 are for composite shells having 4, 8, 16, and 24 axial stiffeners 
with various stiffener heights (1.0", 2.0", 3.0"). The plot of Figure 34 is a comparison 
of the buckling load for a shell having 24 axial stiffeners as a function of stiffener height 
for the various theories used in this study. The results indicate that the discrete LVVTR 
yields more conservative (lower) buckling results than the smeared approaches. More- 
over, as the stiffener height increases the difference between the discrete and the smeared 
approaches increase. In addition, as the number of stiffeners increases the difference 
between the discrete and smeared approaches increases. This could be due to some lo- 
calized stiffener buckling which occurs in the discrete stiffener analyses, but cannot be 
accounted for when using a smeared approach. At this time this localized buckling can- 
not be predicted directly by the discrete method. 

Results for quasi-isotropic [ —45/45/90/0]$ shells with internal (eccentric) ring stifTeners 
subjected to external pressure are presented in Tables 11 and 12 and in Figure 35. The 
results in Tables 6 and 7 are for composite shells having 5, 10, and 25 ring stiffeners 
subjected to external pressure for various stiffener heights (0.5", 1.0", 1.5"). The plot in 
Figure 35 is for buckling pressure versus stiffener height for a cylindrical composite shell 
having 25 ring stiffeners. As can be seen by the results in Tables 1 1 and 12 and in Figure 
35 the discrete LWTR yields more conservative buckling results than the smeared ap- 
proach. The difference becomes more pronounced as the number of stiffeners increases 
and as the stiffener height increases. As mentioned previously this is probably due in 
part to some localized stiffener buckling that is accounted for only in the discrete anal- 
ysis. 
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Table 9. Analysis of [ -45/45/90/0]j 0.2' Thick Circular Cylindrical Shell with Axial Stiffeners Sub- 
jected to Axial Compression (R-85", L- lOtT) - Jones Sroeared/LWTR Discrete. 


No. of Stiffanan 
Stiffener Hefr* 

(mam) 

/ Smeared' 
(feaufacb) 

LMflR 

Oiacreta 

(hi/Mi) 

Error 

Analysis 

Unriltm* 

1799 

3M5, Hm 1 

1793 

M-1 3, N-1 

•0.3 3H 

lwtto 

4/(i.<rj 

2192 

M-1.N-10 

2191 

M-1. N-10 

-q.5 OH 

LWTTO 

• /( 1.0*) 

2255 

ftMI.N-10 

2221 

M-1. N-10 

•1.5391 

LWTTO 

1C / (1.0*) 

2 SCI 

l<M. M-10 

2334 

M-1, M.10 

•2.01 M 

LWTTO 

24/ (1.0*) 

I 

250C 

M-1. M-10 

244C 

M-1. M-10 

•2.4591 

LWTTO 

4/(2^*) 

2423 

M-1.N-10 

2423 

M-1.N-9 

a ok 

lwtto 

• /( 2.0*) 

27 1C 

hUI. N-10 

2543 

CM. M-10 

-6.3791 

LWTTO 

1C / (2.0*) 

329C 

M-1, N-10 

3010 

M-1.N-10 

•9.50H 

LWTTO 

24 / (2.0*) 

3CC7 

M-1, N-10 

3411 

M-1.N-11 

•13*3691 

LWTTO 

4/(3.tr> 

2*40 

CM.N-10 

270* 
M-1. M-1 1 

*4.1 CK 

LWTTO 

• / (3.0* ) 

3743 

M-1, N-10 

317* 

CM.N-11 

-17.7CK 

LWTTO 

ic/(3.«r) 

532* 

M-1. N-10 

423* 

CM.N-11 

•2S.72M 

LWTTO 

24 / (3.0*) 

CM2 

M-1, N-10 

5213 

M-1.M-11 

•31.8691 

LWTTO 


M - numbar of uU halfwaves 
N - numbar of eircumfarantial halfwaves 

UWno - Layer-wise Laminate Theory of Rod* with Oaaatt Stiffmn 
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Table 10. Analysis of [—45/45/90/0]$ 0.2* Thick Circular Cylindrical Shell with Axial Stiffeners 
Subjected to Axial Compression (R*85*, L* 100*) - Reddy Smeared/LNVTR Discrete. 


Mo. of Stiffeners / 
Stiffener Haight 
<***■) 

Smaara d/ 
Rad&Appreadi 

LWTR 

(bayindi) 

Error 

Lowest 

AnWyee 

Utstiffened 

1793 

RMS, N-1 

1783 

ftfelS. N-1 

ao% 

Stmt 

4 /(I .Of) 

2177 

ftfel.N.10 

2181 

ta-i.N-io 

♦a 18% 

Reddy 

• /0.<T) 

2247 

ftfel.lfelO 

2221 

SMI, N-10 

-1.17% 

LWTRD 

1C / (1.0*) 

2W 

RM,N-10 

2334 

ftfel. IfelO 

•2*31% 

LWTRD 

24/ (1.0*) 

2521 

M-1.IM0 

2440 

Mil, N-10 

•3*24% 

LWTRD 

4/ (2.0*) 

2422 

M-1, N-10 

2423 
ftfel, ftfet 

♦ 0 . 04 % 

Reddy 

8 / (2.0*) 

2734 

M-1. N-10 

2543 

RM.N-10 

•7.51% 

LWTRD 

1«/(2-0f) 

3354 

M-1, N-10 

3010 

M, 18*10 

•11.43% 

LWTRD 

24/ (2.0*) 

IH4 

M-1, N-10 

3411 

M-1, N-1 1 

•18*21% 

LWTRD 

4/(S.0f) 

2858 

ftfel.lfelO 

27C8 

ftfel. N-1 1 

-8.75% 

LWTRD 

8/(s.cr ) 

8787 

M-1, N-10 

3178 

M-1, N-1 1 

•19.48% 

LWTRD 

18/ (3.0*) 

5458 

14-1 v N-10 

4233 
M-1, N-1 1 

•28.74% 

LWTRD 

24 / (3.<T) 

7001 

M-1. N-10 

5219 

M-1, N-1 1 

•35.88% 

LWTRD 


M - number of axial halfwave* 

N - nurrtMr of dr cum far an till halfwave* 

LWTRD • L a yer -wi*e Lamina ta Theory of Reddy with Discrete S ti ffener s 
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Table 1 1. Analysis of [ — < 45/45/90/0]$ 0.2* Thick Circular Cylindrical Shell with Ring Stiffeners Sub- 
jected to Lateral Pressure (R “ 85', L ” 100") - Jones Smeared/LWTR Discrete. 


No. of Stiffeners / 
Stiffener Height 
(inches) 

Smeared/ 
Jones Approach 
Bucking Press. 
(D*> 

(LWTR) 
Discrete 
Buddha Press. 

(DSi) 

Error 

Lowest 

Analysis 

Unstiffened 

1.63 

M-1.N-11 

1.61 

M-1.N-11 

-1.1% 

LWTRD 

5/(0.5") 

5.57 

M-l.N-9 

5.00 

M-l.N-9 

-11.4% 

LWTRD 

10/(0.5") 

8.62 

XU1.N-8 

7.49 

M-l.N-9 

-153% 

LWTRD 

2S/(0.5*) 

1533 
M»1, N-8 

12.68 
M-1, N-8 

-223% 

LWTRD 

5/(1.0“) 

18.14 
M-1, N-7 

14.02 

M-1, N-8 

•29.4% 

LWTRD 

10/(1 -0") 

28.91 

M-1.N-7 

21.42 
M-1, N-6 

-35.0% 

LWTRD 

25/(1 -0") 

5138 

*M,N-6 

3639 

M-1.N-5 

•41.6% 

LWTRD 

5/0.5") 

3936 
M-1, N-6 

28.46 

RU1.N-7 

•38.3% 

LWTRD 

1 

1 

10/0-5") 

6236 
M-1, N-6 

28.46 
M-1, N-7 

-44.9% 

LWTRD 

25/(1 .5") 

108.19 
M-1, N-€ 

70.00 

KU1.N-6 

-54.6% 

LWTRD 


M - number of axial halfwaves 
N • number of circumferential halfwaves 

LWTRD • Layer-wise Laminate Theory of Reddy with Discrete Stiffeners 
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Table 12. Analysis of [ -45/45/90/0], 0.2" Thick Circular Cylindrical Shell with Ring Stiffeners Sub- 
jected to Lateral Pressure (R-85', L- 100") - Reddy Smeared/LWTR Discrete. 


|No. of Stiffanera 
jSttffsnsr Haight 
1 (Inches) 

/ Smeared/ 
Reddy Approach 
Bucking Prase. 

m 

(LWTR) 
Discrete 
Bucking Press. 

fosi) 

Error 

Lowest 

Analysts 

1 Uhsttffsned 

1.61 

hM,N-ll 

1.61 

M-1.N.11 

0.0% 

LWTRD 

I 5/(0J") 

050 

M-1.N-9 

540 

H-1.N-9 

• 29 . 9 % 

LWTRD 

I 10/(05-) 

10.03 
M»1, N.0 

7.49 

M-1.N-9 

-34.0% 

LWTRD 

25/(05") 

18.10 

*M,N-7 

12.68 

M-1.N-8 

-42.7% 

LWTRD 

I 5/(14-) 

19.85 | 

M-1.N-7 

1442 

M-1.N-8 

-41.6% 

LWTRD 

10/(14") 

32.00 

M-1.N.7 

21.42 
M-1, N-6 

-49.4% 

LWTRD 

25/(1 4") 

55.71 

kUl.N-6 


•53.5% 

LWTRD 

5/(1.S") 

4143 

M-1.N-6 

28.46 

hUI.N-7 

-47.0% 

LWTRD 

1Q/(1.S") 

66.35 

M-l.N-6 

28.46 
M-l. N-7 

-54.9% 

LWTRD 

25/(1 .5") 

1 13-27 
M-l, Na6 

7040 

hU1.N-6 

•61.8% 

LWTRD 


M - number of udat halfwaves 
N m number of tircumfarendal halfwave s 

LWTRD - Layer-wise Laminate Theory of Roddy with Discrete Stiffeners 
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Analyses of quasi-isotropic [ -45/45/90/0]* shells with eccentric, internal stifTeners was 
conducted next to compare the discrete LWTR with a finite element solution. The finite 
element program CSM Testbed [187,188] was used to analyze the geodesic cylindrical 
shells. Nine-node assumed-natural-coordinate strain (ANS) C° (transverse-shear 
deformable) shell elements were used to model both the shell and the stifTeners. In these 
buckling analyses and in the finite element analyses to follow, geodesically stiffened 
shells having 1x12, 1x16, 2x12, and 2x16 unit cells as described in Figure 36 are studied. 
The finite element model will be made of a unit cell of one of the aforementioned 
geodesic cylinders. A typical finite element model of a unit cell is shown in Figure 37. 
The cylinders were subjected to axial compression and stiffener heights of 0.5", 1.0", 
1.5.", and 2.0" were used for these analyses. The results of this study are presented in 
Tables 13-16. The 2x12 geodesic shell model yields the closest agreement between the 
two analytical methods and the results for this model are plotted in Figure 38. The 
LWTR discrete results and the Testbed finite element results show good agreement. The 
LWTR discrete method yields more conservative buckling loads than the finite element 
method except for the 2x16 shell and the 0.5" and 2.0" stifTener heights in the 1x12 shell. 
The maximum difference in the buckling loads is for the 1x16 shell where the LWTR 
discrete results are 9.3% more conservative for the 1.0" and 1.5" stiffener heights than 
the finite element method. The difference between the LWTR discrete method and the 
finite element method can be attributed to the fact that the LWTR method neglects the 
out-of-plane and the torsional stiffnesses of the stiffeners. Also, the shell is comprised 
of a quasi-isotropic laminate and although small, the orthotropic stiffnesses, 
Ci6, Ci 6 , C 3 6, C 45 , are present. The LWTR discrete approach assumes these values are 
zero and this may result in a slight change in the buckling load. Nevertheless, a good 
correlation of the discrete and finite element buckling results does exist. 
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Glxl 2 


Glxl 6 



Figure 36. Geodesically stiffened shell configurations. 
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I i 8 ure 37. Finite dement mesh used for the geodesie buckling sn.ljsi. (unit cell ni.h 20.20 mesh). 


Table 13. Analysis of [ -45/45/90/0]j 0.2' Thick Circular Cylindrical Shell with Geodesic Stiffeners 
Subjected to Axial Compression (R“85', L“ 100" , 1x12 Geodesic Shell Model). 


Stiffener Height 
(inches) 

Testbed FEA 
(Ibs./inch) 

LWTR 

Discrete 

(lbs./inch) 

Error 

Lowest 

Analysis 

Unstiffened 

1759 

1793 

M-15, N-1 

+1.9% 

FEA 

0.5 

2105 

2119 

+0.7% 

FEA 

1.0 

2172 

2139 

-1.5% 

LWTRD 

1.5 

2232 

2211 

-0.9% 

LWTRD 

2.0 

2284 

2382 

+4.3% 

FEA 


M - 1 , N- 1 0 for all LWTRD Results 

1x12 Geodesic Model 
a - 23.99 Degrees 

M - number of axial halfwaves 
N - number of circumferential halfwaves 

LWTRD - Layerwise Laminate Theory of Reddy with Discrete Stiffeners 
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Table 14. Analysis of [ -45/45/90/0]j 0.2' Thick Circular Cylindrical Shell with Geodesic Stiffeners 
Subjected to Axial Compression (R-85', L« 100', 1x16 Geodesic Shell Model). 


Stiffener Height 
(inches) 

Testbed FEA 
(lbs./inch) 

LWTR 

Discrete 

(lbs./inch) 

Error 

Lowest 

Analysis 

Unstiffened 

1759 

1793 

M-15, N-1 

+1.9% 

FEA 

0.5 

2225 

2121 

-4.3% 

LWTRD 

1.0 

2368 

2148 

-9.3% 

LWTRD 

1.5 

2488 

2257 

-9.3% 

LWTRD 

2.0 

2595 

2494 

-3.9% 



LWTRD 


M - 1 , N - 1 0 for all LWTRD Results 
1x16 Geodesic Model 
a- 18.46 Degrees 


M - number of axial halfwaves 
N - number of circumferential halfwaves 

LWTRD - Layerwise Laminate Theory of Reddy with Discrete Stiffeners 
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Table 15. Analysis of [ -45/45/90/0]s 0.2' Thick Circular Cylindrical Shell with Geodesic Stiffeners 
Subjected to Axial Compression (R“85", L - 100" ,2x12 Geodesic Shell Model). 


Stiffener Height 
(inches) 

Testbed FEA 
(Ibs./inch) 

LWTR 

Discrete 

(lbs./inch) 

Error 

Lowest 

Analysis 

Unstiffened 

1759 

1793 

M-15, N-1 

+1.9% 

FEA 

0.5 

2156 

2129 

-1.2% 

LWTRD 

1.0 

2193 

2146 

-1.8% 

LWTRD 

1.5 

2250 

- 

2189 

-2.7% 

LWTRD 

2.0 

2289 

2286 

-0.1% 

1 

LWTRD 


M - 1 , N - 1 0 for all LWTRD Results 
2x12 Geodesic Model 
a- 41.67 Degrees 


M - number of axial halfwaves 
N - number of circumferential halfwaves 

LWTRD - Layerwise Laminate Theory of Reddy with Discrete Stiffeners 
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Table 16. Analysis of [ -45/45/90/0]* 0.2' Thick Circular Cylindrical Shell with Geodesic Stiflcners 
Subjected to AxiaJ Compression (R-85', L- 100', 2x16 Geodesic Shell Model). 


Stiffener 

Height 

(inches) 

Testbed FEA 
(Ibs./inch) 

LWTR 

Discrete 

(Ibs./inch) 

Error 

Lowest 

Analysis 

Unstiffened 

1759 

1793 

M-15, N-1 

+1.9% 

FEA 

0.5 

2049 

2130 

+4.0% 

FEA 

1.0 

2078 

2146 

+3.3% 

FEA 

1.5 

2122 

2189 

+3.2% 

FEA 

2.0 

2192 

2286 

+4.3% 

FEA 


M ■ 1 , N ■ 1 0 for all LWTRD Results 
2x1 6 Geodesic Model 
a = 33.73 Degrees 


M • number of axial halfwaves 
N - number of circumferential halfwaves 

LWTRD - Layerwise Laminate Theory of Reddy with Discrete Stiffeners 
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Stiffener Height (inches) 



Analysis of cross-ply [0/90/90/0] circular cylindrical shells with eccentric axial, ring, and 
geodesic stiffeners was also performed. The same material properties, shell thickness, 
stiffener thickness, and geometries as those of the quasi-isotropic case were used. 

Results for the cross-ply shells with external, eccentric axial stifTeners subjected to axial 
compression are shown in Tables 17 and 18 and in Figure 39. Once again the discrete 
approach yields more conservative buckling results than the smeared approach. As the 
stiffener height increases the difference between the LWTR discrete and the smeared 
approaches becomes larger as expected. The results plotted in Figure 40 are for 24 axial 
stifTeners at various stiffener heights for the smeared and discrete approaches. 

The results for internally ring stiffened cross-ply cylinders subjected to external pressure 
are provided in Tables 19 and 20 and in Figure 40. The buckling pressure predicted by 
the LWTR discrete approach is much lower than that predicted by the smeared ap- 
proaches. The difference is more pronounced as the stiffener height increases. The plot 
shown in Figure 40 is for a cylinder having 25 internal ring stiffeners. 

Analyses of cross-ply shells with internal geodesic stiffeners was performed to compare 
the LWTR discrete and the CSM Testbed finite element results. The same models and 
geometries used in the quasi-isotropic analysis were also used for the geodesic cross-ply 
analyses (see Figures 36 and 37). The results are for axial compression and are presented 
in Tables 21-24. The 1x12 geodesic shell model yields the closest agreement between the 
two solutions and a plot of the buckling results for this shell is shown in Figure 41. The 
maximum difference in the buckling loads is for the 1x16 shell where the LWTR discrete 
results are 13.4% more conservative for the 1.5' stiffener height than the finite element 
results. As seen from Tables 21 and 23 the results begin to diverge at the 2.0" stiffener 
heights. The difference between the LWTR discrete method and the finite element 
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Table 17. Analysis of [0/90/90/0] 0.2' Thick Circular Cylindrical Shell with Axial Stiffeners Subjected 
to Axial Compression (R “ 85', L - 100') - Jones Smeared/LWTR Discrete. 


No. of Stiffeners / 
Stiffener Hw^it 
(Inches) 

Smeared' 

(fet/tet) 

LOTI 

(toL/indO 

Eirer 

Lowest 

AnWyWs 

Unstiffened 

1033 

M-4.N-17 

1030 

l»M4,N-17 

•a.3% 

LOTTO 

4/<1.0r) 

1379 

M-1. N-10 

1391 

tMI.IMO 

♦0.4% 

Jones 

• /O-ff’) 

1424 

»M.N-10 

ISM 

M-1, N-10 

•left* 

LOTTO 


1513 

M>1, N-10 

1475 

M-1, N-10 

•2.7% 

lwtro 

24/<1.<T) 

1805 

IM,N-10 

1550 

M-1. N-10 

•3.5% 

LOTTO 

4 /ilXT) 

1530 

H«1 f M»10 

1511 

M-1.N-11 

•1-3% 

LOTTO 

9/ (2.(7) 

1726 

14-1. N-10 

1954 

M-1, N-10 

•4-4% 

LOTTO 

19/ (2.0*) 

2109 

hM.N-10 

1972 

M-1. N-10 

•£9% 

LWTRO 

24 / (2.CT) 

2491 

M-1. N-10 

2299 

RM.N-11 

•9-5% 

LOTTO 

4/ (2.0*) 

2049 

M-1, N-10 

1755 

M>1. N-10 

•15-3% 

LOTTO 

9/ (3.0*) 

2745 

M-1 f N-10 

2197 

M-1, N-10 

•25.7% 

LWTRO 

16 / (3 xn 

4120 

RM.N-10 

2579 

hM.N-11 

•38.3% 

LOTTO 

24/ (3.0*) 

54SO 

M-1.N-10 

3729 

l^1«N-11 

-46-2% 

LWTRO 


M - numbar of atdal halfwaves 
N - number of droumfarantial halfwaves 

LWTRO - Layer-wise Laminate Theory of Reddy with Ossete Stiffeners 
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Table 18. Analysis of [0/90/90/0] 0.2- Thick Circular Cylindrical Shell with Axial Stiffeners Subjected 
to Axial Compression (R-85", L- 100") - Reddy Smeared/LWTR Discrete. 


No. of Stiffenen / 
Sdf 

(Indus) 

Reddy Approach 

(— l/Mi) 

LET— 

(ta*Ancft) 

Eiror 

Lowest 

AfWyW 

Unstiffened 

1030 

M-4, N-17 

1030 

M-4.N-17 

ao% 

LWTRD 

4/(i^r) 

1374 

*M,N-10 

1311 

RM t N-10 

♦05% 

Reddy 

• /(I XT) 

1427 

3M. N-10 

ISM 

M>10 

•2.0% 

LWTRD 

1i/(1.0») 

1530 

RM. N-IO 

- 

1475 

Mai* N-10 

•3.7% 

LWTRD 

24/ (1.0*) 

1131 

M-1, N-10 

1550 

M-1, N-IO 

•5J% 

LWTRD 

4 / (2.<T) 

1523 

bM.N-10 

1511 

bM.N-11 

«4et% 

LWTRD 

• /(2.0T) 

1S4f 

M»1 9 N-10 

1154 

RM,N-10 

-11. IN 

LWTRD 

n/( 2 .<r) 

2354 

RM.N-10 

1072 

M-1. N-10 

•11.5% 

LWTRD 

24/(2.0T) 

2153 

M»1, N-10 

2211 

M»1,N-11 

•1SJR 

LWTRD 

4/(8.0P) 

2071 

M-1, N-10 

1755 

bM.N-10 

•107% 

LWTRD 

• /(3.C) 

2107 

M-1 f H*10 

2117 

RU1.N-10 

•21.3% 

LWTRD 

11/ (3.0") 

4240 

RM1.N-10 

2575 
►4-1. M- 11 

-42.4% 

LWTRD 

24/(3.CT) 

5127 

RU1, N-10 

3720 

M-1.N-11 

•501% 

LWTRD 


M • number of axial halftravas 
N - number of circumferential halteavos 

UffTEO • Layer wise Laminate Theory of Reddy wtm OMarete Stiffeners 
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Figure 39. Axial buckling for a 24 axial stiffener shell model ([0/90/90/0] layup; R m 85*, L *= 100*; 1.0* x 0.2* stiffeners). 



Table 19. Analysis of [0/90/90/0] 0.2* Thick Circular Cylindrical Shell with Ring Stiffeners Subjected 
to Lateral Pressure (R - 85', L - 100*) - Jones Smear ed/LWTR Discrete. 


No. of Stiffeners 
Stiffener Height 
(inches) 

/ Smeared/ 
Jones Approach 

Budding Press. 
(o2) 

(LWTR) 
Olscrete 
Bucking Press. 

(PSD 

Error 

Lowest 

Analysis 

Unstiffened 

0.89 

*M,N-12 

0.88 

M-1.N-12 

-0.9% 

LWTRD 

5/(05") 

4.67 

M-1.N-9 

4.08 

M-1.N-9 

-14.4% 

LWTRD 

10/(0.5") 

72 1 

M-1, N-8 

6^1 

M-1.N-9 

-16.0% 

LWTRD 

25/(0.5") 

12 27 
M-1, N-7 

10.63 
M-1, N-7 

-21.0% 

LWTRD 

5/(1^-) 

15.36 
M-1. N-7 

11.73 
M-1, N-7 

-30.9% 

LWTRD 

10/(1.0") 

23.83 

hU1.N-6 

18^5 

M-1.N-6 

-30.6% 

LWTRD 

25/(1 «0") 

42.49 

M-1.N-5 

29.69 

kM,N-6 

-43.1% 

LWTRO 

5/(1 .S”) 

32.73 

M-1.N-6 

23.414 

M-1.N-6 

-39.8% 

LWTRD 

10/OJ-) 

49.92 
M-1, N-S 

36.19 

M-1.N-6 

-37.9% 

LWTRD 

25/(1 .5") 

9072 

M-1.N-4 

55.75 

kM.N-S 

•62.7% 

LWTRD 


M • number of ajdal halfwaves 
N - number of tircumfarential halfwaves 

IWTRD - Layer-wise Laminate Theory of Reddy with Discrete Stiffeners 
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Table 20. Analysis of [0/90/90/0] 0.2' Thick Circular Cylindrical Shell with Ring Stiffeners Subjected 
to Lateral Pressure (R“85', L* 1 00") - Reddy Smeared/LWTR Discrete. 


No. of Stiff anen 
Stiffener Hei^tt 
(Inches) 

Smeared/ 
Reddy Approach 
Budding Press. 

(os!) 

(LWTR) 
Dricrete 
Budding Press. 

(Dif) 

Error 

Lowest 

Analysis 

Unstiffened 

0.88 

M-1,N-12 

048 

M-1.N-12 

0.0* 

LWTRD 

5/(04*) 

542 

M-1, N-8 

448 

M-1.N-9 

-35.1 * 

LWTRD 

10/(0.5") 

847 

M-1, N-8 

641 

M-1, N-8 

-38.0* 

LWTRD 

25/(0.5”) 

15.15 

RM.N-7 

10.63 

M-1.N-7 

•424* 

LWTRD 

S/(UD 

1743 

M-1.N-7 

11.73 

M-1. N-8 

•45.1* 

LWTRD 

10/(1. 0") 

26.02 

M-1.N-6 

1845 

RM.N-6 

-42.6* 

LWTRD 

25/(14-) 

4545 
M-1 N-5 

29.69 

M-1.N-6 

-52.7* 

LWTRD 

5/(1. 5«) 

35.17 

M-1.N-6 

23.41 

bUI.N-6 

•504* 

LWTRD 

10/(14-) 

5240 
M-1, N-5 

36.19 

M-1.N-6 

•45.9* 

LWTRD 

25/(14") 

93.76 

M-1.N-4 

55.75 

RM.N-S 

•68.2* 

LWTRD 


M - number of arid halfwaves 
N • number of dreumferential halfwaves 

LWTRD - Layer-wise Laminate Theory of Reddy with Discrete Stiffeners 


Results 


15S 




150 



Results 


1S6 


Figure 40. Buckling pressure for . 25 ring stiffener shell model ([0/90/90/0] l.yup; R - g S «, I. - |00~; |.r * 0.2' stiffeners). 



method can be attributed to the fact that the LWTR method neglects the out-of-plane 
and the torsional stifTnesses of the stiffeners. The difference starts to become more 
pronounced as the stiffener heights increase. This is due to the fact that as the stiffener 
height increases the effects of the out-of-plane and the torsional stiffnesses on the global 
buckling results increase. Except for the 1x16 geodesic shell model there appears to be 
a good correlation of the buckling results, especially for the 1.0" and 1.5" stiffener 
heights. 

Finally, a study of the buckling of geodesically stiffened cylinders subjected to external 
pressure was conducted. A 1x12 geodesic shell model was selected which has the same 
geometry and dimensions as used for the axial compression analyses. A comparison 
between the LWTR discrete method and the finite element method (CSM Testbed) was 
made for both the cross-ply case, [0/90/90/0], and the quasi-isotropic case, 
[ -45/45/90/0]s. The results for the cross-ply [0/90/90/0] shell is found in Table 25. 
The LWTR results correlate fairly well with the finite element solutions especially for the 
unstiffened case and for lower stiffener heights. The quasi-isotropic results are presented 
in Table 26. The results indicate that at lower stiffener heights the buckling pressures 
for both analytical methods are reasonably close, but as the stiffener height increases the 
buckling pressures tend to diverge. Neglecting the orthotropic stiffnesses, 
Ci 6 , C 26 , Cm, C 45 , definitely must have a major impact on the stiffness of the layerwise 
quasi-isotropic laminates when the shells are sbjected to external pressures. Another 
difference between the LWTR discrete method and the finite element method can be 
attributed to the fact that the LWTR method neglects the out-of-plane and the torsional 
stiffnesses of the stiffeners. The difference becomes larger as the stiffener height in- 
creases because as the stiffener height increases the effects of the out-of-plane and the 
torsional stiffnesses on the global buckling results become more prominent. 
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The typical CPU time for the Testbed finite element buckling analyses is 620 seconds. 
The CPU time for the layerwise discrete method for one buckling mode (m, n) is 410 
seconds. Consequently, the CPU times for the layerwise discrete method can become 
large, perhaps 8- TO hours or more, if a sweep of a large number of buckling modes (m, 
n) is made in order to determine the minimum eigenvalue. 


5.2 LWTR/Testbed Finite Element Stress Analysis Comparison 

A stress analysis comparison of the LWTR and CSM Testbed finite element codes was 
made for geodesically stiffened shells. A 1x12 geodesically stiffened shell as shown in 
Figure 36 with 1.0" high by 0.2" thick internal stiffeners was selected as the comparison 
model. The shell geometry consists of a radius of 85" and length of 100". In lieu of 
modeling the entire cylinder, symmetry conditions were employed and an analysis of a 
unit cell was made, see Figures 42 and 43. The loadings were employed via the appli- 
cation of uniformly applied end displacements. Three laminate layups were studied in 
this analysis: [0/90/0]; [45/ - 45/45/ - 45]; and [60/ - 60/0/ - 60/60]. The material 
properties used are given in Table 7. The ply thickness used in these analyses is 0.100". 
Analyses were performed with 0° orthotropic stiffeners. Applied end displacements of 
0.01 were used to generate the compressive loads. 

The CSM Testbed elements described in section 5.1 were used for this analysis. The fi- 
nite element model shown in Figure 42 describes the finite element mesh and the 
boundary conditions used for this analysis. The Testbed finite element model uses 256 
nine-node assumed-natural-coordinate strain (ANS) shear deformable shell elements to 
model the the shell and 32 nine-node plate elements to model the stiffeners. 
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Table 21. Analysis of [0/90/90/0] 0.2" Thick Circular Cylindrical Shell with Geodesic Stiffeners Sub- 
jected to Axial Compression <R“85", L m 100", 1x12 Geodesic Shell Model). 


Stiffener Height 
(inches) 

T estbed FEA 
(Ibs./inch) 

LWTR 

Discrete 

(Ibs./inch) 

Error 

Lowest 

Analysis 

Unstiffened 

1053 

1030 

M-15, N-1 

-2.2% 

LWTRD 

0.5 

1347 

1326 

-1.6% 

LWTRD 

1.0 

1570 

1549 

-1.3% 

LWTRD 

1.5 

1664 

1702 

+2.2% 

FEA 

2.0 

1743 

1953 

+10.8% 

FEA 


M - 1 , N- 1 0 for all LWTRD Results 

1x12 Geodesic Model 
a - 23.99 Degrees 

M - number of axial halfwaves 
N - number of circumferential halfwaves 

LWTRD - Layerwise Laminate Theory of Reddy with Discrete Stiffeners 
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Table 22. 


Analysis or [0/90/90/0] 0.2' Thick Circular Cylindrical Shell with Geodesic Stiffeners 
Jected to Axial Compression (K - 85', L - 100', 1x16 Geodesic Shell Model). 


Sub- 


Stiffener Height 
(inches) 

Testbed FEA 
(Ibs./inch) 

LWTR 

Discrete 

(lbs./inch) 

Error 

Lowest 

Analysis 

Unstiffened 

1053 

1030 

M-15, N-1 

-2.2% 

LWTRD 

0.5 

1357 

1327 

-2.2% 

LWTRD 

1.0 

1768 

1537 

-13.1% 

LWTRD 

1.5 

1980 

1715 

-13.4% 

LWTRD 

2.0 

2150 

2016 

-6.2% 

LWTRD 


M - 1 , N- 1 0 for all LWTRD Results 
1x16 Geodesic Model 
a - 1 8.46 Degrees 


M - number of axial halfwaves 
N - number of circumferential halfwaves 

LWTRD - Layerwise Laminate Theory of Reddy with Discrete Stiffeners 
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Table 23. Analysis of [0/90/90/0] 0.2' Thick Circular Cylindrical Shell with Geodesic Stiffeners Sub- 
jected to Axial Compression (R“85', L* 100", 2x12 Geodesic Shell Model). 


Stiffener Height 
(inches) 

Testbed FEA 
(Ibs./inch) 

LWTR 

Discrete 

(Ibs./inch) 

Error 

Lowest 

Analysis 

Unstiffened 

1053 

1030 

M-15, N-1 

-2.2% 

LWTRD 

0.5 

1505 

1588 

+5.2% 

FEA 

1.0 

1724 

1764 

+2.3% 

FEA 

1.5 

1802 

1908 

+5.6% 

FEA 

2.0 

1838 

2052 

+10.4% 

FEA 


M « 1 , N- 1 0 for all LWTRD Results 
2x1 2 Geodesic Model 
a - 41.67 Degrees 

M - number of axial halfwaves 
N • number of circumferential halfwaves 

LWTRD - Layerwise Laminate Theory of Reddy with Discrete Stiffeners 
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Table 24. Analysis of [0/V0/90/0] 0.2* Thick Circular Cylindrical Shell with Geodesic Stiffeners Sub- 
jected to Axial Compression (R-85*, L- 100*, 2x16 Geodesic Shell Model). 


Stiffener Height 
(inches) 

Testbed FEA 
(Ibs./inch) 

LWTR 

Discrete 

(lbs./inch) 

Error 

Lowest 

Analysis 

Unstiffened 

1053 

1030 

M-15, N-1 

-2.2% 

LWTRD 

0.5 

1658 

1574 

-5.2% 

LWTRD 

1.0 

1937 

1763 

-9.0% 

LWTRD 

1.5 

2052 

1957 

-4.6% 

LWTRD 

2.0 

2180 

2206 

+1.2% 

FEA 


M - 1 , N- 1 0 for all LWTRD Results 
2x1 6 Geodesic Model 
o « 33.73 Degrees 

M - number of axial halfwaves 
N « number of circumferential halfwaves 

LWTRD - Layerwise Laminate Theory of Reddy with Discrete Stiffeners 
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Figure 41. Axial buckling results for a 1x12 geodesic shell model ([0/90/90/0] layup; R - 85', L - 100*; I.O” x 0.2' stiffeners). 



Table 25. Analysis of [0/90/90/0] 0.2' Thick Circular Cylindrical Shell with Geodesic Stiffeners Sub- 
jected to Lateral Pressure (R«85', L- 100', 1x12 Geodesic Shell Model). 


Stiffener 

Height 

(inches) 

Testbed FEA 

(psi) 

LWTR 

Discrete 

(psi) 

Error 

Lowest 

Analysis 

Unstiffened 

0.89 

0.88 

M-1, N-11 

-1.1% 

LWTRD 

0.5 

1.22 

1.17 

-4.0% 

LWTRD 

1.0 

1.88 

1.67 

-11.2% 

LWTRD 

1.5 

2.40 

2.17 

-9.6% 

LWTRD 

2.0 

2.65 

2.33 

-12.1% 

LWTRD 


M - 1 , N - 11 for all LWTRD Results 
1x12 Geodesic Model 
a *23.99 Degrees 


M - number of axial halfwaves 
N * number of circumferential halfwaves 

LWTRD - Layerwise Laminate Theory of Reddy with Discrete Stiffeners 
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Table 26. Analysis of [ -45/45/90/0]* 0.2' Thick Circular Cylindrical Shell with Geodesic StilTeners 
Subjected to Lateral Pressure (R-85', L- 100', 1x12 Geodesic Shell Model). 


Stiffener- 

Height 

(inches) 

Testbed FEA 
(psi) 

LWTR 

Discrete 

(psi) 

Error 

Lowest 

Analysis 

Unstiffened 

1.63 

1.61 

M-1, N-11 

-1.2% 

LWTRD 

0.5 

2.09 

1.90 

-9.1% 

LWTRD 

1.0 

2.79 

2.09 

-25.1% 

LWTRD 

1.5 

3.80 

2.20 

-42.1% 

LWTRD 


M - 1 , N - 1 1 for all LWTRD Results 
1x12 Geodesic Model 
a=23.99 Degrees 

M ■ number of axial halfwaves 
N = number of circumferential halfwaves 

LWTRD - Layerwise Laminate Theory of Reddy with Discrete Stiffeners 
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The LWTR finite element model described in Figure 43 shows the finite element mesh 
and the boundary conditions used for this analysis. A more refined mesh is used in 
critical regions such as the stiffener intersection. The model employs 256 four-node 

LWTR shell elements to model the shell and 32 two-node LWTR beam elements to 
model the stiffeners. 

The first analysis considered here is for a [0/90/0] shell with orthotropic (0°) stiffeners. 
The plot in Figure 44 shows the axial stress, <r„, for all layers at x = ^ plotted along 
the nondimensional circumference of the unit cell, The stresses calculated from the 
two analyses are in good agreement away from the point of the stifiener intersection. 
Near the stiffener intersection i.e. ± , the LWTR stresses are less compressive 

than the stresses calculated by the Testbed analysis. For the bottom 0° layer the axial 
compressive stress is 3.2% (120 psi) less compressive than the Testbed axial stress. This 
small difference is within reason and some of the difference can be attributed to the fact 
that a different type of element is used for the stiffeners in each of the models and a small 
difference in the behavior of the stiffener intersection intersection is being observed. The 
axial stresses m the 90° layer as shown in Figure 44 are in good agreement except at the 
location of the stiffener intersection where the LWTR stress, <?„, is 10.6% (30 psi) less 
compressive than the Testbed axial stress at that location. The LWTR axial stress at the 
stiffener intersection in the top 0° layer is 3.0% (104 psi) less compressive than the 
Testbed axial stress. It is apparent from this analysis that the stiffener intersection is 
slightly more compliant (less stiff) for the LWTR method when compared with the re- 
sults generated from the Testbed analysis. 

The second analysis considered here is that of an angle ply laminate, 
[45/ - 45/45/ - 45], with orthotropic stiffeners. The axial stresses for the individual 
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layers are presented in Figures 45 and 46. The axial compressive stresses in the individual 
layers are in fairly close agreement. The largest difference in stresses occurs at the 
stiffener intersection location for the top — 45° layer shown in Figure 46 where the 
LWTR are 10.5%' (40 psi) less compressive than those calculated from the Testbed 
method. This difference in stresses at the stiffener intersection can be attributed to the 
difference in the stiffener intersection stiffness of the two analytical methods. 

The final analysis is that of a symmetric 5 layer quasi-isotropic shell laminate, 
[60/ — 60/0/ — 60/60], with orthotropic stiffeners. The results of this analysis are pre- 
sented in Figures 47-49. There is good agreement of the axial stresses for all layers. The 
axial stresses in the 60° and —60° layers are 5-6% (20-25 psi) smaller in the LWTR along 
the circumference of the shell. The largest percentage difference in compressive stress 
occurs at the location of the stiffener intersection in the top 60° layer where the LWTR 
axial stress is 10.4% (35 psi) smaller than the Testbed stress, see Figure 49. In the 0° 
layer the Testbed results are an average of 2.5% (80-100 psi) smaller than the LWTR 
results. Thus, for the LWTR analysis the 0° layer is carrying slightly more compressive 
load and the 60° and —60° layers do not carry quite as much compressive load when 
compared with the Testbed results. Overall there is good agreement of the axial stresses 
for all layers. 

The stress analysis comparison here was conducted to help verify the stress analysis ca- 
pabilities of the LWTR finite element program. Several lamination schemes were con- 
sidered to accomplish this task. A good correlation exists between the LWTR analyses 
and the Testbed analyses. Small differences in stresses do occur at the siffener inter- 
section. These differences are not major and are due to the difference in stiffener inter- 
section response measured by the two analysis methods. 
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The CPU time for the Testbed finite element stress analysis is 125 seconds. Run times 
for the layerwise theory vary depending upon the number of shell laminate layers, num- 
ber of nodes in the mesh, and the type of elements used (4 or 9 nodes). A three layered 
shell, [0/90/0], model (289 nodes) with one layered beams has 3567 active degrees of 
freedom. A mesh of 256 four node layerwise shell elements and 32 two node layerwise 
beam elements will have a half-bandwidth of 285. The CPU times for this model are 139 
seconds for a linear analysis and 730 seconds for a geometrically nonlinear analysis (4 
iterations to converge). If the mesh is changed to 64 nine node layerwise shell elements 
and 16 three node layerwise beam elements with a half-bandwidth of 555, the CPU times 
increase to 736 seconds for a linear analysis and 4251 seconds for a nonlinear analysis 
(5 iterations to converge). A six layered shell, [45/90/0]* with 289 nodes using one 
layered beams has 6168 active degrees of freedom. A mesh of 256 four node layerwise 
shell elements and 32 two node layerwise beam elements will have a half-bandwidth of 
456. The CPU times for this model are 752 seconds for a linear analysis and 3261 sec- 
onds for a geometrically nonlinear analysis (4 iterations to converge). Increasing the 
bandwidth has a large influence on the CPU time necessary to run the layerwise finite 
element analyses. Furthermore, layerwise elements are not practical elements to use if 
a postbuckling analysis must be conducted due to the large run times. 


5.3 Displacements and Interlaminar Stresses in Geodesically Stiffened Shells 


In this study it was desired to conduct analyses of geodesically stiffened shells using the 
LWTR in order to determine the trends that have major effects on the transverse dis- 
placements and the interlaminar stresses. Variables such as the laminate layup and 
thickness, stiffener height, stiffener orientation, stiffener angle, cell geometry, cell length, 
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Figure 42. Unit cell finite element mesh anil boundary conditions for the Testbed stress analysis (256 elements). 
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Hgure 43. Unit cell finite element mob and boundary conditions for the I.W’IK stress analysis (256 elements). 



figure 44. LWTR and Testbed 1x12 geodesic shell axial stresse for [0/90/0] laminate: a) bottom 0° pi); b) 90° ply; c) top 0° ply ( 









b) bottom *45* ply 

45 -T-a jss^s, stirs ;? t , *•& ■— 
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b) top -45° ply 


Figure 46. LWTR and Tcstlicd 1x12 geodesic shell axial stresses for top layers of 
[45/ - 45/45/ - 45] laminate: a) top 45° ply; b) top -45° ply (x - L/2). 
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Figure 48. 


L\VTR and Testbed 1x12 geodesic shell 
[60/ - 60/0/ - 60/60] laminate (x - L/2\ 


stresses for 0° layer of 
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Fijture 49. LWTR and Testbed 1x12 teodesic *holl .„;.i .. 

[60/ - 60/0/ - 60/60] laminate: a) top 60* ply; b) top -60* plyoT" L/2). * 0P ****” ° f 
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were varied when performing this study. Concentration here is placed upon the action 
of the shells under pure compressive loads and combined loads generated with the addi- 
tion of internal pressure. In addition, a geometrically nonlinear analysis was performed 
to determine the effects on the displacement and stress fields. A base line design using 
a 1x12 geodesically stiffened shell with a nominal radius of 85", shell length of 100", 1.0" 
high by 0.2" thick orthotropic internal stiffeners, and a [0/90/0] shell laminate was used 
for a large number of the analyses. A simple shell laminate, [0/90/0], was used to per- 
form many of the comparison studies in order to keep the number of degrees of freedom 
manageable and in particular to keep the bandwidth of the global stiffness matrix from 
becoming excessively large. The bandwidth can become extremely large when analyzing 
a large number of nodes and laminate layers using the layerwise theory. The LWTR fi- 
nite element model used in these analyses was described previously in Figure 36. The 
model consists of 289 nodes, 256 four-node layerwise shell elements, and 32 two-node 
beam elements. 


5.4 Displacement Field in Geodesically Stiffened Shells 

A study of the transverse displacement field along the circumference of the unit cell at 
x « — was performed for various shell parameters. The transverse displacements are 
the largest and most interesting displacements for the structures being analyzed in this 
research. The shells are subjected to compressive loading via the application of applied 
end displacements of 0.01" and to a combined loading consisting of applied end dis- 
placements of 0.01" and an internal pressure of 10 psi. The transverse displacements 
presented in this study are nondimensionalized using the following expression: 
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X 100 


(5.4) 


_ vvZT. h 2 
w = ■■■ __ — 

A(,RL 

Mere 5? is the nondimensional transverse displacement, w is the transverse displacement, 
h is the total laminate thickness, £, is the modulus in the 1 direction, N x is the applied 
load, R is the shell radius, and L is the shell length. The first study considered here is 
the efiect of changing the laminate layup when the shell is subjected to compressive 
loading. The nondimensional displacement results presented in Figure 50 are for six 
di/Terent laminate layups: [0/90/0], [45/ - 45/45/ - 45], [60/ - 60/0/ - 60/60], 

[45/90/0]$, [ -45/45/90/0]$, and [30/ - 30/0/ - 30/30]. The [0/90/0], [ -45/45/90/0]$, 
and [30/ - 30/0/ - 30/30] laminates show the largest variations in nondimensional 
transverse displacement 234%, 155%, and 113% respectively from the edge of the unit 
cell to the stilTener intersection at = 0 .5. The [45/ - 45/45/ - 45], 

[60/ - 60/0/ - 60/60], [45/90/0]$ yield respective changes in nondimensional transverse 
displacement along the circumference of the shell of 45%, 52%, 100%. The 

[45/ - 45/45/ - 45] and [30/ - 30/0 / - 30/30] laminates yield the maximum nondimen- 
sional transverse displacement along the circumference of 1.484 and 1.305 respectively 
for the given geometry and loading conditions. Thus, in order to avoid the largest 
transverse displacements and variations in transverse displacements in geodesically stiff- 
ened composite shells it is best to avoid designs containing cross-ply laminates, angle 
ply laminates, and laminates containing +30° or —30° plys. 

The next study concentrates on the effects of changing the geodesically stiffened shell 
geometry and the stiffener orientation angle on the transverse displacement Held. A 
[0/90/0] 1x12 geodesically stiffened shell is used in this analysis. The shells are subjected 
to uniform compressive loading via applied end displacements of O.Or. The stiffener 
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Figure SO. Nondimensional transverse displacements for a 1x12 geodesically stiffened shell as a function 
of the laminate stacking sequence under compressive loading (R ■ 85*, L " 1 00*). 
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orientation angle, a, relative to the shell is shown in Figure 1 and Figure 14. Figure 51 
shows the nondimensional transverse displacement as a function of the geodesic shell 
geometry and the subsequent stifTener orientation angle. The variation in the stifTener 
orientation angles were obtained by changing the cell geometry from G 1x12 to Glxl6 
in Figure 51a and by changing the shell length from 50" to 200" in Figure 51b. The 
nondimensional transverse displacement field presented in Figure 51a shows that as the 
cel! geometry is changed from Glxl2 (a = 24°) to Glxl6 (a = 18.5°) the maximum non- 
dimensional transverse displacement decreases by 72%. The maximum transverse dis- 
placement shows the largest decrease, 41%, when the cell geometry is changed from 
Glxl2 (a = 24°) to G 1x1 3 (a = 22.3°). When the cell geometry is changed from Glxl5 
(a « 19.6°) to Glxl6 (a = 18.5°) the change in maximum transverse displacement be- 
comes less significant (2.5%). The results presented in Figure 51b show the transverse 
displacement as a function of the stifTener angle by changing the shell length. The trends 
predicted here indicate that as the shell length is decreased from 200" to 50" the maxi- 
mum transverse displacement decreases by 52% even though the stifTener orientation 
angle is increased from 12.5° to 41.7°. Thus, decreasing the shell length has a much 
greater effect on the shell stiffness than does increasing the stiffener orientation angle. 
The results presented in Figure 51 indicate that the response of geodesic shells is a 
stronger function of the shell geometry i.e. the number of cells around the circumference 
and the cell length and a weaker function of the stiffener orientation angle. 

The effect of increasing the stiffener height on the nondimensional transverse displace- 
ment for both compressive loading and combined loading is shown in Figure 52. A 
[0/90/0] 1x12 geodesically stiffened shell is used in this analysis. Shown in Figure 52a 
is a plot of the nondimensional transverse displacement along the nondimensional shell 
circumference. As expected the nondimensional transverse displacements decrease as a 
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Figure SI. Nondimensional transverse displacements fur geodesically stiffened shells as a function of the stiffener orientation: 
a) changing cell geometry; b) changing cell length. 





function of the stiffener height. As the stiffener height is increased the area of the 
stiffeners is increased and consequently the applied load N x increases because the load 
is a function of the applied displacements. In this study the applied displacement are 
constant and thus- increasing the area increases the applied load. The transverse dis- 
placements are also a strong function of the stiffener bending stiffness and to some ex- 
tent the axial stiffness of the stiffener. Thus, the transverse displacements are a function 
of the applied load, which increases as the stiffener height is increased, and the stiffness 
of the stiffeners. Shown in Figure 52b is the effect of the stiffener height on the nondi- 
mensional transverse displacements along the shell circumference when the shell is sub- 
jected to combined loads. The decrease in the nondimensional transverse displacement 
appears to be fairly uniform. In this case the applied load N x is a function of the applied 
displacements, load area, and the internal pressure. The load area is a function of the 
increasing stiffener height. Therefore, the normalized transverse displacements will be a 
complex function of the stiffener height, stiffener stiffness, and the internal pressure. 

The nondimensional transverse displacements as a function of the shell laminate thick- 
ness for a [0/90/0] 1x12 geodesically stiffened shell are shown in Figure 53. The laminate 
thicknesses studied here are 0.15”, 0.30”, 0.45”, and 0.60”. The results for compressive 
loading, shown in Figure 53b, reveals that the nondimensional transverse displacements 
decrease by 40.6% when increasing the shell thickness from 0.15” to 0.60". The maxi- 
mum normalized transverse displacements for the 0.15" and 0.30” thicknesses are close 
together. This probably occurs because the nondimensional transverse displacement 
used in this study is a function of the shell thickness squared. The actual transverse 
displacements differ by 88%. The maximum difference in the actual transverse dis- 
placements is between the 0.15" and 0.60" shell thicknesses and is 558%. The results for 
the combined loading condition as a function of shell laminate thickness is shown in 
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*) compressor loading; b) combined loading. 





Figure 53b. As in the case for compressive loading, increasing the shell laminate thick- 
ness causes the nondimcnsional transverse displacement to decrease. The maximum 
difference in the nondimcnsional transverse displacement is 40.3% when increasing the 
shell thickness from 0.15" to 0.60". The actual transverse displacement difference be- 
tween the 0.15" and 0.60" thicknesses is 386%. The transverse displacements are deli- 
nitely a function of the shell laminate thickness. 

Shown in Figure 54 are the nondimensional transverse displacement comparisons of 
linear and geometrically nonlinear analyses Tor [0/90/0] 1x12 geodesically stiffened 
shells. The results for the compressive loading load is found in Figure 54b. The nondi- 
mensional transverse displacement results indicate that using a geometrically nonlinear 
analysis does yield a suffer structure and consequently slightly smaller transverse dis- 
placements than the linear analysis. The maximum difference in the nondimensional 
transverse displacements between the linear and nonlinear analysis is 8.9%. The results 
for a combined loading are shown in Figure 54b. The addition of internal pressure yields 
larger differences and variations in the nondimensional transverse displacements between 
the linear and nonlinear analyses than does the case of compressive loading only. As 
can be seen from Figure 54b the distribution and magnitude of the nonlinear results are 
much different than the linear results. The maximum difference in the displacements is 
28.7% for the combined loading condition. Therefore, it can be concluded that the ge- 
ometric nonlinearities become more significant with the addiuon of internal pressure. 

The effect of increasing the shell radius is shown in Figure 55. The nondimensional 
transverse displacements for shell radii of 85". 170", and 255" subjected to an applied end 
compression developed through applied displacements for a [0/90/0] laminate are de- 
scribed in Figure 55a. As can be seen, as the sheU radius is increased the maximum 
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nondimensional transverse displacement is reduced by 52.8% when the shell radius is 
increased from 85" to 170". When the shell radius is increased from 85" to 255" the 
maximum nondimensional transverse displacements are reduced 73%. Figure 55b shows 
the nondimensional transverse displacements when an internal pressure of 10 psi is 
added to the applied end compression loading for varying shell radii. When the radius 
is increased from 85" to 170" the maximum transverse increases 47%. Increasing the 
shell radius from 85" to 255" increases the transverse displacement by 55.5%. The ad- 
dition of internal pressure yields large differences and variations in the nondimensional 
transverse displacements as a function of the shell radius. As can be seen from Figure 
55b, the distribution of the nondimensional transverse displacements yields a maximum 
displacement at the stiffener intersection for the 85" shell, almost a uniform variation 
of the displacements for the 170" shell, and a minimum displacement at the stiffener 
intersection for the 255" shell. A possible explanation for this phenomena is that for a 
constant shell thickness, as the shell radius is increased the ratio of the stiffener stiffness 
to the shell stiffness to resist pressure induced deflections increases. Moreover, as the 
shell radius is increased the load on the shell due to internal pressure increases and the 
stiffeners exert a larger influence on the structural response. This becomes evident with 
the addition of pressure loading. In this case the 85" shell has the smallest stiffness ratio 
and thus with the addition of pressure the stiffeners tend to deflect more. The 170" shell 
has an almost equal ratio and therefore a uniform displacement field is observed. The 
255" shell has the largest stiffness ratio and thus with the addition of pressure the shell 
deflects more away from the stiffener intersection. 

Shown in Figure 56 is a comparison of a [0/90/0] 1x12 geodesically stiffened shell, a 
[0/90/0] 1x12 axial/ring stiffened shell, and a [0/90/0] unstiflened shell. The axial/ring 
stiffened shell internal axial and ring stiffeners consisting of 1.0" orthotropic stiffeners. 
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The axial/ring stiffened finite element model with appropriate boundary conditions is 
shown in Figure 57. The results for compressive end loads are shown in Figure 56a. The 
geodesically stiffened shells show the largest transverse displacements of the shells 
studied here. The nondimcnsional transverse displacements for the geodesically stiffened 
shells are 91% greater than the displacements for axial/ring stiffened shells and arc 90% 
greater than the displacements of unstiffened shells. The results for combined loading 
is shown in Figure 56b. The geodesically stiffened shells exhibit the largest nondimen- 
sional transverse displacements. The geodesically stiffened shells yield nondimensional 
transverse displacements 46.3% larger than the axial/ring stiffened shells and are 37.1% 
larger than the unstiffened shells. The axial/ring stiffened shells show that the shell 
stiffness at the stiffener intersection is much greater than that of the geodesically stiff- 
ened shell when internal pressure is applied. In this case the displacements away from 
the stiffener intersection are larger than those at the intersection in much the same way 
as the displacement field generated via the 255” shell shown in Figure 55b. 

The last displacement field analysis involves studying the effects of combined loading 
on geodesically stiffened shells. Figure 58 shows the results of these loadings on 0.075", 
0.15", and 0.30" [0/90/0] laminates with various stiffener heights. The results show that 
the maximum nondimensional transverse displacement occurs at the stiffener inter- 
section until the stiffeners reach a critical height at which point the stiffeners are suffi- 
ciently stiff in bending to prevent large transverse displacements at the stiffener 
intersection. At this point the shell displacements away from the stiffener intersection 
becomes larger than the displacements at the stiffener intersection. The stiffener height 
at which this occurs is also a function of the laminate thickness. For example, the 
transverse displacement away from the stiffener intersection exceeds the transverse dis- 
placement at the stiffener intersection for a stiffener height of 2.0” for a 0.075" laminate, 
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1'igurc 57. Unit cell finite element mesh and boundary conditions of the LWTR stress analysis of the axial and ring stiffened shell (256 elements). 






40" for a 0.15" laminate, and 6.0" for a 0.30" laminate. Therefore, under the action of 
compressive load and internal pressure the stiffeners act to decrease the transverse dis- 
placements at the stiffener intersection, but only when the stifTeners are sufficiently deep 
compared to the laminate thickness do the shell transverse displacements away from the 
stiffener intersection exceed those at the stiffener intersection. 

This study of the transverse displacement field trends has yielded some interesting re- 
sults. The shell laminate layup will have a major factor on the transverse displacement 
field. Changing the stiffener height, shell geometry, shell radius, and shell laminate 
thickness all have a major impact on the structural response of the stiffened shells. A 
geometric nonlinear analysis does not yield major changes in the displacement field for 
compressive loading, but nonlinearity is significant when pressure is added to the load- 
ing. The geometry of the geodesically stiffened shell has a significant impact on the 
displacement field. When subjected to combined compression and internal pressure the 
transverse displacements of the shell away from the stiffener intersection do not exceed 
those at the stiffener intersection until deep stiffeners, a large radius, or an axial/ring 
stiffened structure are used. One of the nice features in using the LWTR finite element 
code to conduct this design analysis is that it is fairly simple and quick to generate new 
models by changing the cell geometry or the stiffener parameters. This is not true for the 
Testbed finite element code where a more time consuming effort is needed to generate 
models that change shell geometries and/or stiffener heights. 
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5.5 Detailed Stress Study 


The results presented in this section will focus on interlaminar stresses, but a few exam- 
ples of the in-plane stresses a xx and a xy over the region of the entire shell will also be 
presented. Interlaminar stresses in the geodesic shells have never been studied in detail 
before. The interlaminar stresses over the entire shell structure for a few specific cases 
will be presented. This will help to determine the regions of peak interlaminar stresses 
and the nature of the stress distribution over the entire region. The interlaminar stresses 
at the critical regions, probably near the stiffener intersection, will be studied. The 
interlaminar stress distribution through the thickness at the critical regions will be 
studied. The effects of the shell laminate layups, laminate thickness, pressure loading, 
stiffener height, shell radii, cell geometry, and geometric nonlinearity on the interlaminar 
stresses will presented. The base line design used in this study is a [0/90/0] 1x12 
geodesically stiffened shell with a shell radius of 85“ and 1.0“ internal orthotropic 
stiffeners. The base line laminate thickness used here is 0.30“. The loadings considered 
in these analyses are applied compressive end loads generated through applied end dis- 
placements of 0.01 on each edge (x = 0, Lx) and combined compressive loads and 
internal pressure (10 psi). Some of the more interesting interlaminar stress results are 
presented in this work. The stresses are nondimensionalized in this study using the fol- 
lowing expression: 


a — 


ghL 

A^R 


(5.5) 


Here a is the nondimensional stress, a is the generated stress, h is the total laminate 
thickness, N x is the applied load, R is the shell radius, and L is the shell length. 
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5.5.1 In-Plane Stress Study 


The nondimensional in-plane stresses, a xx and a xy , will be discussed in this section. The 
in-plane stress, a xx , of the inner layer for a [0/90/0] 1x12 geodesically stiffened shell with 
1.0" internal stiffeners subjected to an applied end load of 0.01" is shown in Figure 59. 
The peak in-plane stresses, a a , occur at the boundaries x = 0, L* and also at the 
stiffener intersection x - -y- , y - -y\ The in-plane stresses, o„, are the largest of the 
six stresses and thus are likely to be the primary contributing stresses to cause failure. 
By viewing Figure 59 it can be said that failure would most likely occur at the boundaries 
(x = 0, Lx) or at the stiffener intersection. Figure 60 is a plot of over the stiffened 
shell structure for the outer layer of a [0/90/0] 1x12 geodesically stiffened shell with 1.0" 
internal orthotropic stiffeners subjected to a compressive end load. The stress distrib- 
ution, over the surface of the shell in the top layer as shown in Figure 60 yields a 
different stress field shape than that generated in the inner layer shown in Figure 61. 
The stresses in the outer layer peak at the boundary comers and at the stiffener inter- 
section. As can be seen from Figures 59 and 60, the in-plane stresses, <?**. for the inner 
layer are more uniformly distributed, particularly at the boundaries than the in-plane 
stresses for the outer layer. One possible explanation for this phenomena is that the 
stiffeners are attached to the inner layer and this reduces the bending of the inner layer 
and in addition some of the load is carried by the stiffeners. This results in more uniform 
stresses in the inner layer. The influence of the stiffeners on the outer layers are evi- 
denced by the fact that the in-plane stresses are lower at the comers of the boundaries 
and at the stiffener intersection. However, at the center of the shell boundaries the in- 
fluence of the stiffeners is not as pronounced and more bending occurs. This results in 
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the outer layer carrying more load at these locations which produces larger compressive 
stresses. 

The in-plane stress results, a xx% for a combined loading on a [0/90/0] 1x12 geodesically 
stiffened shell are shown in Figures 61 and 62. The addition of 10 psi internal pressure 
to the geodesically stiffened shell produces significant changes in the stress distribution, 
o xx , in the inner layer shown in Figure 61 when compared with the inner layer for 
compressive loading shown in Figure 59. As the internal pressure is increased these 
differences between pure compression and combined loading will become much larger. 
The pressure produces a much wider variation of the in-plane stresses and also changes 
the peak stresses at the edges and the stiffener intersection. The in-plane stress at the 
stiffener intersection for the inner layer subjected to combined loading is reduced by 52% 
from the inner layer stress generated from compressive loading. Figure 62 shows the 
stress distribution, for the outer layer for the base line design. The boundaries at x 
= 0, Lx are fairly stiff having v - w - 0 boundary conditions at these locations and 
thus the addition of internal pressure results in large compressive bending stresses at the 
x boundaries. Away from the boundaries the pressure tends to reduce the compressive 
stresses by as much as 84% at the stiffener interior. Thus, adding pressure has a signif- 
icant influence on the in-plane stresses of the shell laminate. 

A plot of the shear stress, in the inner layer for the base line design subjected to 
combined loading is shown in Figure 63. The shear stresses, <?„, yield a skew-symmetric 
nature with the value of the shear stress being 0 at the stiffener intersection. The results 
for the application of compressive loading yield the same general shape as the results for 

combined loading, but the values of the stresses are about a factor of 10 smaller. Those 
results are not included here. 
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Figure 59. 


Surface plot of for the inner layer of a [0/90/0] 1x12 geodesically stiffened shell under 
compressive loading. 
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Figure 60. 
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Figure 61. Surface plot of <F„ for the inner layer of a [0/90/0] 1x12 geodesically stiffened shell under 
combined loading. 
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Figure 62. 


Surface plot of ?„ for the outer layer of a [0/90/0] 
combined loading. 


1x12 geodesicaily stiffened shell under 
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Figure 63. Surface plot of?,, for the inner layer of a [0/90/0] 1x12 geodcsically stiffened shell under 
combined loading. 
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5.5.2 Interlaminar Normal Stress Study 


The distribution of the interlaminar normal stress, e„. over the shell region for the outer 
layer or the base line design subjected to an applied compressive end load is shown in 
Figure 64. As can be seen, the interlaminar normal stresses peak at the x boundaries 
and also at the stiffener intersection. Along the center line y - h. away from the 

boundaries the interlaminar normal stresses are larger than those stresses over the re- 
mamder of the shell away from the surface. This indicates that the stiffener intersection 
has an influence on the interlaminar normal stresses along the line y - h.. This couW 
be due to the fact that the transverse displacements do peak at the stiffener intersection. 
Adding internal pressure to the preexisting compressive load yields an interlaminar 
normal stress distribution in the outer layer as described in Figure 65. Here the 
interlaminar normal stresses are significantly greater than for the case of end com- 
pression. The general pattern of the stress distribution is the same as that for 
compressive end loading only. The interlaminar normal stresses peak at the stiffener 
intersection and the interlaminar normal stresses being largest along the line y - -L. 
Pressure does have a significant influence on the interlaminar normal stresses by in- 
creasing the interlaminar normal stresses by almost a factor of 4. 

Nondimensional interlaminar normal stresses near the stiffener intersection through the 
thickness of the shell laminate for various geometries and loadings are presented in this 
section. The interlaminar normal stresses generated via combined loading are 400% 
larger than the interlaminar normal stresses generated from compressive loading and 
thus only combined loading conditions are studied in this section. Figure 66 is a plot 
of the interlaminar normal stresses, 3,„ near the stiffener intersection for 1x12 
geodesically stiffened shells having I. O' interna, orthotropic stiffeners for various shell 
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Figure 64. Surface plot of a„ for the outer layer of a [0/90/0] 1x12 geodcsically stiffened shell under 
compressive loading. 
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laminates subjected to combined loading. Figure 66 shows that the laminate stacking 
sequence has a large impact on the interlaminar normal stress. The 
[60/ - 60/0/ — 60/60] and the [—45/45/90/0]$ layups show the largest interlaminar 
normal stresses. For example, the maximum nondimensional interlaminar normal 
stresses in the [60/ — 60/0/ — 60/60] layup is 66% larger than the maximum nondimen- 
sional interlaminar normal stress in the [45/90/0]$. The order of decreasing maximum 
nondimensional interlaminar normal stress stresses in the laminates are 
[30/ - 30/0/ - 30/30], [ -45/45/90/0]$, [45/ - 45/45/ - 45], [60/ - 60/0/ - 60/60], 

[0/90/0], and [45/90/0]$. Hence, laminates such as the [45/90/0]$, [0/90/0], and 
[60/ - 60/0/ — 60/60] are preferable for use in keeping the stresses from becoming 
significantly large. 

The effects of conducting a geometrically nonlinear analysis on the stresses, a„, for a 
combined loading is shown in Figure 67. The nondimensional interlaminar normal 
stresses for the geometrically nonlinear analysis are 47.1% less than the stresses devel- 
oped from the linear analysis. From the displacement field study shown in Figure 54b 
the nonlinear analysis generates smaller displacements than the linear analysis and ob- 
viously this results in lower strains and then subsequently lower stresses. Thus, when 
pressure loading is included on the structure a geometrically nonlinear analysis will yield 
significantly different nondimensional interlaminar normal stresses. The geometrically 
nonlinear analysis softens the structure and this reduces the displacements and subse- 
quent stresses. 

The effect of increasing the stiffener height on the nondimensional interlaminar normal 
stress for combined loading is shown in Figure 68. The nondimensional interlaminar 
normal stresses show a uniform decrease as the stiffener height is increased. The maxi- 
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Figure 66. ThrougMhe-thickness distribution of 7 a for Glxl2 shell near the stiffener intersection for 
various shell laminates under combined loading. 
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mum differences in a zz between the 1" and 2", 3", 4", and 5" stiffeners are 6.6%, 18.4%, 
27 . 9 %, and 37.5% respectively. The interlaminar normal stress decreases as the stiffener 
height increases as expected because the bending stiffness of the stiffeners increases as 
a function of the cube of the stilfener height. This increase in bending stiffness tends to 
decrease the transverse displacements and consequently the interlaminar normal stress. 
It is interesting to note that the height of the stiffener does not affect the overall shape 
of the stress through the thickness. 

The efTect of changing the cell geometry on the nondimensional interlaminar normal 
stress for [0/90/0] geodesically stiffened shells subjected to combined loading is shown 
in Figure 69. Geodesic cell geometries of 1x10, 1x12, 1x14, and 1x16 are considered in 
this analysis. As can be seen from Figure 69 and the displacement results shown in 
Figure 51, increasing the number of cells around the circumference causes the shells to 
become significantly stiffer because the number of stiffeners in the structure is increased. 
This lowers both the transverse displacements and the interlaminar normal stresses. For 
example, as the cell geometry is increased from 1x10 to 1x12, 1x14, and 1x16 the non- 
dimensional interlaminar normal stresses are reduced 35%, 51%, and 85% respectively. 

Shown m Figure 70 are the nondimensional interlaminar normal stress results for in- 
creasing the shell radius for [0/90/0] shells subjected to combined loading. The shell 
radii considered here are 85", 170”, 255". The maximum difference in nondimensional 
interlaminar normal stress is 56% when increasing the shell radius 2 times from 85" to 
170 and the difference is 88% when increasing the shell radius 3 times from 85” to 255". 
As can be seen there are differences in the nondimensional interlaminar normal stress 
distribution. From the results shown in Figure 70 and the transverse displacement re- 
sults of Figure 55b it is observed that as the shell radius is increased the influence of the 
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Through-the-thickness distribution of a„ for Glxl2 shell near the stiffener intersection for 
varying stiffener heights under combined loading. 
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Figure 69. Through-the-thickness distribution or* tt for Glx!2 shell near the stiffener intersection for 
vorying the ceil geometry under combined loading. 





stiffener on the displacement and interlaminar normal stress fields increases. This is es- 
pecially evident when pressure is applied to the shell structure. Moreover, as the shell 
and stiffener radius are increased and the shell thickness remains constant, the ratio of 
the stiffener stiffness to the shell stiffness increases and the stiffeners' influence on the 
response of the shell is increased. 

The nondimensional interlaminar normal stresses for [0/90/0] 1x12 geodesically stiffened 
shells subjected to combined loading for variations in the shell laminate thickness are 
shown in Figure 71. As expected, changes in shell thickness have a direct effect upon 
the interlaminar normal stresses. The nondimensionalized interlaminar stresses are a 
strong function of the shell laminate thickness. The maximum differences in nondimen- 
sional interlaminar normal stresses between the 0.15' and the 0.30', 0.45”, and 0.60' 
laminates are 78%, 89%, and 97% respectively. The bending stiffness of the shell is a 
function of the laminate thickness cubed. The smaller the laminate thickness the more 
the shell will deflect under pressure loading resulting in larger stresses. 

The nondimensional interlaminar normal stresses for a 1x12 geodesically stiffened shell, 
1x12 axial/ring stiffened shell, and an unstiffened shell are shown in Figure 72. The re- 
sults presented here are for [0/90/0] shells with 1.0' internal orthotropic stiffeners sub- 
jected to combined loading. As can be seen the nondimensional interlaminar normal 
stresses generated in the geodesically stiffened shell is 67.8% larger than those generated 
in the axial/ ring stiffened shell and 81.7% larger than those generated from the unstiff- 
ened shell. From Figures 56b and 72 it becomes apparent that the geodesic stiffeners 
tend to push the stiffener intersection outward which results in larger transverse dis- 
placements and stresses near the stiffener intersection than those generated by the 
axial/ring stiffened shell system or the unstiffened shells. Apparently, the geodesic 
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Figure 71. Through-the-thicknes* distribution of for Glx12 shell nemr the stiffener intersection for 
increasing shell laminate thickness under combined loading. 
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stiffeners at the stiffener intersection are constrained so that the u and v displacements 
at that location are 0 because of symmetry and symmetric loading and boundary condi- 
tions. However, because of this constraint the stiffeners do exhibit large displacements 
in the transverse direction due to the compliance of the stiffener intersection. In addi- 
tion, the distribution of the interlaminar normal stresses for the axial/ring stiffened shell 
and the geodesically stiffened shell are different. While both stiffened shells exhibit peak 
stresses in the bottom (inner) layers of the laminate where the stiffeners are attached, the 
axial/ring stiffened shell does not exhibit another peak in the interlaminar normal stress 
in the top (outer) layers as does the geodesically stiffened shells. Therefore, it can be 
concluded that because the geodesically stiffened shells produce significantly more dis- 
placements at the stiffener intersection than the axial/ring stiffened results as seen in 
Figure 56b this results in larger interlaminar normal stresses through the entire shell 
laminate for the geodesically stiffened shells. The axial/ring stiffened shells are very stiff 
at the stiffener intersection which results in smaller displacements at the stifiener inter- 
section than away from the intersection (see Figure 56b). Therefore, the influence of the 
axial/ring stiffened structure is to cause peak interlaminar normal stresses in the inner 
layers of the shell, but because the transverse displacements away from the stiffener 
intersection are larger than at the stiffener joint the influence of the axial/ring stiffeners 
does not extend to the outer layers of the shell. 


The interlaminar normal stresses are influenced by the laminate stacking sequence, ge- 
ometric nonlinearity, stifiener height, cell geometry, shell radius, shell laminate thickness, 
and the type of shell structure (geodesic or axial/ring stiffened). The shape and magni- 
tude of the nondimensional interlaminar normal stress is definitely influenced by the 
laminated stacking sequence. Changing the shell geometry, shell laminate thickness, shell 
radius, and conducting a geometrically nonlinear analysis all have an impact in the 
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structural response and the interlaminar normal stresses. Comparing a 1x12 geodesically 
stiffened shell, 1x12 axial/ring stiffened shell, and an unstiffened shell shows that the 
interlaminar normal stresses generated in a geodesically stiffened shell are much larger 
than those generated in other types of shells. However, the interlaminar normal stresses 
are still an order of magnitude less than the in-plane stresses. The addition of internal 
pressure increases the interlaminar normal stresses and therefore large increases in pres- 
sure may cause the interlaminar normal stresses to contribute to the failure of the 
geodesically stiffened shells. 


5.5.3 Interlaminar Shear Stress Study 

The interlaminar shear stresses, are the interlaminar shear stresses having the largest 
magnitudes for geodesically stiffened shells and will be studied in this work. The 
interlaminar shear stress distribution, <?„, for the outer layer base line design subjected 
to a compressive loading is shown in Figure 73. The interlaminar shear stress surface 
plot yields a skew- symmetric stress distribution. The shear stress do not peak at the 
stiffener intersection, but rather peak about 3.5” from the stiffener intersection. The 
interlaminar shear stress, over the shell for a combined applied compressive load and 
an internal pressure of 10 psi is shown in Figure 74. This reveals that the interlaminar 
shear stresses yield similar type of behavior and peak at the same location as shown in 
Figure 73. The difference in nondimensional interlaminar shear stress between the 
compressive and combined loading is only 28.6%. Thus, a large porti * of the 
interlaminar shear stresses, a xl , are generated by the compressive rather tha e pres- 
sure loads. The combined loading case does produce larger nondimensional interlaminar 
shear stresses by 40% at the x boundaries, x = 0, L,. Combined loading will be used 
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to study the interlaminar shear stresses in order to be consistent with the interlaminar 
normal stress analyses. 

Shown in Figure 7.5 are the nondimensional interlaminar shear stresses, a zl , at the lo- 
cation of the peak stresses as indicated from the surface plots of Figures 73 and 74 for 
various shell laminates. The nondimensional interlaminar shear stresses for the 
geodesically stiffened shells shown in Figure 75 are developed via combined loading. 
Clearly due to the influence of the stiffeners and the stacking sequence there is no dis- 
tinct pattern for the shear stresses a xr The [30/ - 30/0/ - 30/30], [0/90/0], and the 
[45/ — 45/45/ — 45] laminates yield the maximum values of the stresses a xt . The 
[45/90/0]$ laminated stiffened shell yields the smallest nondimensional interlaminar 
shear stresses. The difference in the nondimensional shear stresses between the 
[30/ — 30/0/ — 30/30] shell and the [45/90/0]$ shell is 84.5%. Obviously the shell lami- 
nate has a definite influence on the interlaminar shear stresses. Laminates such as the 
[ —45/45/90/0]$, [45/90/0]$, and [60/ — 60/0/ — 60/60] are preferable for use in keeping 
the a xt stresses from becoming significantly large. 

In Figure 76, the interlaminar shear stresses, a xt , are compared for a linear and a ge- 
ometrically nonlinear analysis subjected to combined loading of the base line design. 
The results show that the nondimensional interlaminar shear stresses produced from the 
linear analysis are 80% larger than those generated from the geometrically nonlinear 
analysis. It is clear that a nonlinear analysis does tend to soften the structure, especially 
when pressure is applied, thus reducing the resulting displacements, strains, and stresses. 

Next, the effect of changing the cell geometry upon the nondimensional interlaminar 
shear stresses for a [0/90/0] laminated shell subjected to combined loading is shown in 
Figure 77. The results presented here are for 1x10, 1x12, 1x14, and 1x16 geodesically 
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Figure 73. 


Surface plot of for the outer layer of a [0/90/01 1x12 
compressive loading. 
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Figure 74. Surface plot of for the outer layer of a [0/90/0] 1x12 geodesically stiffened shell under 
combined loading. 
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Figure 75. Throufh-the-thJcliijes. distribution of c„ for Clxl2 shell at the critical region for various 
shell laminates under combined loading. 
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G1x12, 0/90/0, R=85", haO.30*. SH=1.0", pi=10 psi, del=0.01* 

G1x12, 0/90/0, R=85\ h=0.30*. SH=1.0', pi=10 psi, defcO.OT (nonlinear) 



Figure 76. Through-the-thickness distribution of for Glxl 2 shell at the critical region for linear and 
geometrically nonlinear analyses under combined loading. 
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stiffened shell models. Obviously as the number of cells around the circumference from 
1x10 to 1x16 the shell becomes stiffer and thus lower interlaminar stresses result. The 
resulting differences in the maximum nondimensional interlaminar shear stresses be- 
tween the 1x10 and 1x12, 1x14, and 1x16 models are 57%, 86%, and 93% respectively. 

Changing the shell laminate thickness for the base line design under combined loading 
has a definite impact on the interlaminar stresses a„ as shown in Figure 78. The results 
presented in Figure 78 are for 0.15", 0.30', 0.45'. and 0.60' shell laminate thicknesses. 
The nondimensional shear stress results for the 0.15" laminate shell thickness are 87% 
greater than those for the 0.30" laminate, 96.2% greater than those for the 0.45' lam.- 
nate, and 98.2% greater than those for the 0.60' laminate. The bending stiffness of the 
shell is a function of the laminate thickness cubed. Under pressure loading the smaller 

laminate thicknesses will definitely deflect more and therefore larger interlaminar shear 
stresses are developed. 

In Figure 79 the effects of increasing the shell radius on the nondimensional interlaminar 
shear stress a xt under combined loading is described. Shell radii of 85", 170", and 255" 
are considered in this analysis. The shell laminate thickness remains a constant 0.30", 
The results presented in Figure 79 clearly show that increasing the shell radius while 
holding the shell laminate thickness constant definitely has an impact on the 
interlaminar shear stresses. The addition of internal pressure is the loading which brings 
out the large variation in the interlaminar shear stress a xt . As discussed in Section 5.4 
concerning Figure 55b, increasing the shell radius while holding the laminate thickness 
constant increases the ratio of the stiffener stiffness to the shell stiffness to resist pressure 
loading. Thus, near the stiffener intersection for the 255" shell the geodesic stiffeners 
have a large influence on the displacement field when pressure is applied and therefore 
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Figure 77. Through-the-thickness distribution of <F„ for GlxI2 shell it the critical region for changing 
cell geometry under combined loading. 
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shcli laminate thickness under combined loading. 
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larger shear stresses result in the bottom laminate layers near the stiffener intersection. 
In fact, the shear stresses change sign through the laminate thickness for the 255" shell. 
There appears to be a balance between the increased stiffness created when increasing 
the shell radius and the increased loading generated via internal pressure by increasing 
the shell radius. Thus, when the shell radius is increased from 85" to 170" the effect of 
the increased stiffness is not overcome by the increased loading from the internal pres- 
sure and thus the interlaminar stresses decrease. However, when the shell radius is in- 
creased from 85" to 255" the increase in loading from internal pressure exceeds the 
increase in shell stiffness and thus the interlaminar shear stresses exhibit a large variation 
through the laminate thickness. 

The last interlaminar shear stress analysis considered here involves 1x12 geodesically 
stiffened shell, 1x12 axial/ ring stiffened shell, and an unstiffened shell as shown in Figure 
80. The results presented here are for [0/90/0] shells with 1.0" internal orthotropic shells 
under combined loading. The results shown in Figure 80 clearly show that the influence 
of the stiffeners on the interlaminar shear stresses when compared with an analysis of 
an unstiffened shell. The maximum difference between the geodesically stiffened shell 
and the unstiffened shell is 93.3% while the maximum difference between the axial/ring 
stiffened shell and the unstiffened shell is 95.5%. Also, the distribution of the 
interlaminar shear stresses for the axial/ring stiffened shell and the geodesically stiffened 
shell are different. While both stiffened shells exhibit peak stresses in the bottom (inner) 
layers of the laminate where the stiffeners are attached, the axial/ring stiffened shell does 
not exhibit another peak in the interlaminar shear stress in the top (outer) layers as does 
the geodesically stiffened shells. Therefore, it can be concluded that because the 
geodesically stiffened shells produce significantly more displacements at the stiffener 
intersection than the axial/ring stiffened results as seen in Figure 56b this results in larger 
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interlaminar shear stresses through the entire shell laminate for the geodesically stiffened 
shells. The axial/ring stiffened shells are very stifT at the stiffener intersection which re- 
sults in smaller displacements at the stiffener intersection than away from the inter- 
section (see Figure 56b). Therefore, the influence of the axial/ ring stiffened structure is 
to cause peak interlaminar shear stresses in the inner layers of the shell, but because the 
transverse displacements away from the stiffener intersection are larger than at the 
stiffener joint the influence of the axial/ring stiffeners does not extend to the outer layers 
of the shell. 

The impact of laminate stacking sequence, geometric nonlinearity, cell geometry, shell 
laminate thickness, shell radius, and shell type (unstifTened, axial/ring stiffened, or 
geodesically stiffened) all influence the interlaminar shear stress. Changing any or many 
of these parameters will result in significant changes in the interlaminar shear stresses. 
However, these shear stresses are in many cases 2 orders of magnitude less than the in- 
plane stresses and these shear stresses' impact on the structural integrity is debatable. 

5.6 First-Ply Failure Analysis 

The purpose of this study is to determine the primary failure trends for geodesically 
stiffened shells. This work is not intended to be a detailed failure study similar to the 
research presented by Reddy and Pandey [159], but rather it should be viewed as a pre- 
liminary study of the failure of geodesically stiffened shells. The laminate strength values 
used in this analysis are given in Table 7. The Tsai-Wu failure criterion discussed in 
section 2.6 will be used to determine shell laminate material failure. Shown in Table 27 
are a comparison of the results for linear material failure and buckling of [0/90/0] 1x12 
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geodesically stiffened shells with 1.0" internal orthotropic stiffeners. Three shell laminate 
thicknesses, 0.15", 0.30", and 1.5", were considered in this analysis. Results were ob- 
tained for pure compression and for shells subjected to combined loading with internal 
pressures of 10 psi and 25 psi. The end displacements were increased incrementally until 
failure occurred. The end loads, N X| were calculated at first-ply failure. The results of 
Table 27 clearly show that for pure compression and combined loading with small to 
moderate internal pressures the shell will buckle before material failure occurs. It is also 
evident that as the pressure is increased from 0 psi to 25 psi the difference between the 
failure and buckling loads decreases rapidly. The buckling load increases as the pressure 
is increased due to the biaxial loading. Because the unit load in the circumferential 
(hoop) direction, N y , is tensile when internal pressure is included, the buckling load, N x , 
increases due to this biaxial loading condition. The failure load decreases as the pressure 
is increased and thus it is concluded that for larger pressures material failure will occur. 

The next study shown in Table 28 is for the analysis of [0/90/0] unstifTened shells, 
[0/90/0] 1x12 geodesically stiffened shells, and [45/90/0]$ 1x12 geodesically stiffened 
shells. A constant shell thickness of 0.30" is used. For these analyses the end displace- 
ment is held constant and the internal pressure is increased until material failure occurs. 
Applied end displacements of 0.0", 0.10", and 0.25" were used as the constant end dis- 
placements. As can be seen material failure occurs at high pressures, (160 psi - 225 psi), 
and this depends upon the shell type and the applied displacements. The failure of the 
geodesically stiffened shells occurs at higher pressures than for the unstifTened shells 
primarily because of the stress concentrations at the stiffener intersection. As discussed 
in sections 5.2, 5.4, and 5.5 the maximum displacements, in-plane stresses, and 
interlaminar stresses occur at the stiffener intersection. Shown in Figure 8 1 is the finite 
element model and the location of first-ply failure marked with a circled X. The corners 
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of the shell and the stiffener intersection are the locations of first-ply failure for the 
geodesically stiffened shells as expected. However, some of the peak stresses at the 
corners may be artificially induced by the boundary conditions. The location of the 
first-ply failure depends upon the laminate layup and the applied end displacements. 
The [0/90/0] shells tend to fail at the stiffener intersection first and the [45/90/0]$ shells 
tend to fail at the shell corners first. The differences in the failure pressures between the 
unstiffened and geodesically stiffened shells range from a minimum for 4.4% for the 
[0/90/0] geodesically stiffened shell with 0.0" end displacements to a maximum of 15.6% 
for the [45/90/0]$ geodesically stiffened shell with 0.0" end displacements. The maxi- 
mum difference in failure pressures between the [0/90/0] and the [45/90/0]$ geodesically 
stiffened shells occurs for applied end displacements of 0.0". The failure pressures of the 
geodesically stiffened shells are the same for applied end displacements of 0.25". In all 
cases there is a large variation in the applied end load. It can be concluded that material 
failure of geodesically stiffened shells will occur for large internal pressures in the vicinity 
of the stiffener intersection or the corners of the shell section. 
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Table 27. Comparison of First- Ply Failure and Buckling Loads for Geodeskally Stiffened [0/90/0] Shells. 
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1x12 Geodesically Stiffened Shells with 1.0" Internal Orthotropic Stiffeners 



Table 2*. First Ply Failure Results for Geodesically Stiffened 
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Figure 81. Location of the first-ply failure in the laycrwise finite element model. 


Chapter 6 

Conclusions and Recommendations 


6.1 Summary and Conclusions 

The purpose of this research was to analyze geodesically stiffened shells using a layerwise 
approach. The literature review, theoretical developments, verifications of the analytical 
method and computer codes, and the analysis of the stiffened shells were all vital to the 

completion of this research. A summary of the major accomplishments of this work 
follow: 

• Extensive literature review including shell theories, buckling of stiffened shells, finite 

element analysis of stiffened shells, failure theories, and calculation of interlaminar 
stresses. 

• Theoretical and computational development of a layerwise discrete Ritz buckling 
procedure. 
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• Buckling analysis of geodesic, axial, and ring stilTened shells with smeared and finite 
element comparisons. 

• Theoretical and computational development and subsequent verification of the 
layerwise shell and stiffener finite elements with applications to geodesically stiffened 
shells and interlaminar stresses. 

The literature review provided the groundwork for this research and any subsequent 
work w'hich naturally follows from this research. All of the analytical developments were 
derived because of a lack of study or knowledge about a particular area or because of 
an interest in expanding the database about a certain topic. 

The layerwise discrete Ritz analytical buckling procedure was developed to prove the 
validity of the layerwise theory for use in the analysis of geodesically stiffened composite 
shells. Attachment of the discrete stiffeners was implemented by using the Lagrange 
multiplier technique. The layerwise analytical buckling results compare well with the 
smeared buckling results and the Testbed finite element results. The layerwise discrete 
analytical method yielded more conservative buckling results than the smeared results 
and the differences in results ranged from 0% to 70%. This is reasonable because dis- 
crete methods should yield more conservative buckling results than the smeared ap- 
proaches. A comparison of the Testbed finite element buckling results with the layerwise 
discrete results shows that in general for quasi-iso tropic and cross-ply laminates the 
layerwise theory yields more conservative results. The difference between the LWTR 
discrete method and the finite element method can be attributed to the fact that the 
LWTR method neglects the out-of-plane and the torsional stiffnesses of the stiffeners. 
The difference starts to become more pronounced as the stiffener heights increase. This 
is due to the fact that as the stiffener height increases the effects of the out-of-plane and 
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the torsional the stiffnesses on the global buckling results increase. The layerwise dis- 
crete method may be employed for axial, ring, or geodesically stiffened shells. This 
method is only directly applicable to specially orthotropic shells. Reasonable buckling 
results should be expected. The layerwise discrete buckling method provides design en- 
gineers with an optional tool in the design of stiffened shells based upon buckling. One 
of the drawbacks of this method is that it is restricted to shells which have boundary- 
conditions that are analytically tractable. Run times for this method can be extreme if 
a sweep of the buckling modes must be made in order to determine the minimum 
eigenvalue. Also, this method is only directly applicable to specially orthotropic shells 

where C l6 = C 26 = C J6 = C* 5 = 0. However, this method has provided good results 
for certain quasi-isotropic materials subjected to in-plane loads. 


The layerwise finite element method for geodesically stiffened shells was developed pri- 
marily to study the displacement and stress fields in geodesically stiffened shells. Of 
particular interest were the interlaminar stresses. Both the layerwise shell and beam el- 
ements were developed for this research. The out-of-plane stiffness of the layerwise 
beam elements was included by using the ratio of the out-of-plane moment of inertia to 
the m-plane moment of inertia. Neglecting the out-of-plane beam stiffness has the 
greatest impact when angle ply laminates are used in the shell or beams. For these cases 
the finite element method does not yield good results unless the out-of-plane beam 
stiffness is included. Developing the beam elements in a layerwise fashion permits the 
beam element degrees of freedom to be assembled directly into the global stiffness ma- 
trix. Thus, no additional constraint equations are necessary. The layerwise finite ele- 
ment program was verified using 10 classical example problems. An additional 
comparison of the layerwise finite element method with the Testbed finite element 
method was conducted for geodesically stiffened shells and several shell lamination 
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schemes. A good correlation of stress results exists between the layerwise analyses and 
the Testbed analyses, less than 10%. The small differences occur at the stiffener inter- 
section and are due to the different stiffener interaction responses measured by the two 
analyses. 

A study of the displacements and the stresses for various geodesically stiffened shells was 
made using the layerwise finite element method. Varying the shell laminate layup, lam- 
inate thickness, stiffener height, stiffener orientation angle, cell geometry, shell radii, and 
shell length were all considered . All changes have an impact on the structural response, 
some more than others. The displacement field is most affected by changing the shell 
laminate layup and the cell geometry. Adding internal pressure to the shell has a major 
influence on the displacement response when compared with pure compression. This is 
most evident for the geometrically nonlinear case, when the shell radius is increased, and 
when axial and ring stiffened shells are considered. When subjected to combined com- 
pression and internal pressure the transverse displacements of the shell away from the 
stiffener intersection do not exceed those at the stiffener intersection until deep stiffeners, 
a large radius, or an axial/ring stiffened structure are used. 

Interlaminar normal and shear stresses for geodesically stiffened shells are only signif- 
icant when internal pressure is added. The interlaminar normal stresses are influenced 
by the laminate layup, geometric nonlinearity, stiffener height, cell geometry, shell ra- 
dius, shell laminate thickness, and the type of shell structure. The shape and magnitude 
of interlaminar normal stress distribution through the thickness is influenced most by the 
laminate layup, shell laminate thickness, and shell type (unstiffened, axial and ring stiff- 
ened, or geodesically stiffened). However, the interlaminar normal stresses are an order 
of magnitude less than the in-plane compressive stresses. The interlaminar shear streses 
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are greatly influenced by the laminate layup, shell laminate thickness, and shell type. 
The interlaminar shear stresses are often two orders of magnitude less than the in-plane 
stresses. The influence of the interlaminar stresses on the structural integrity of stiffened 
shells is small for pure compression and combined loading when the internal pressure is 
small. 


The failure analysis reveals that unless internal pressure is applied the geodesically stiff- 
ened shells will buckle before they experience material failure. Increasing the internal 
pressure can create a failure scenario. For large internal pressures failure of geodesically 
stiffened shells initiates at the stiffener intersection where the largest displacements and 
stresses occur. 

The layerwise finite element method provides a useful analytical tool to study the struc- 
tural response of geodesically stiffened shells. The layerwise method eliminates the finite 
element aspect ratio problem of traditional 3-D finite elements. Also, the layerwise finite 
element code was written so that a preprocessor is not needed to generate a large model 
and thus it is easier to change variables such as the shell radius, shell length, cell geom- 
etry, and stiffener height than for the Testbed models. However, for large problems the 
size of the bandwidth can hamper the solution by causing excessive run times and not 
converging to the correct solution. 

6.2 Recommendations 


The recommendations include expanding the existing analytical tools and augmenting 
the analyses. The layerwise Ritz method should be expanded to include calculation of 
the displacements, strains, and stresses for simply supported cylindrical shells. Including 
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the out-of-plane and torsional stiffnesses of the stiffeners used in the Ritz method can 
also be included in the future. Employment of layerwise beam elements in the Ritz 
method may also be useful. Also, perhaps the layerwise discrete approach could be ex- 
tended to other tractable boundary conditions and linear vibration analyses. Develop- 
ment of robust postprocessors for both the analytical and finite element codes will 
provide the user with plots of the deformed shapes, eigenvectors, and stress contours. 
Including the torsional stiffnesses of the layerwise beam elements by assuming a dis- 
placement distribution through the thickness of the beam could improve finite element 
method. Improvement and/or additional finite element equation solvers could help cir- 
cumvent the bandwidth problems for large models. Parallel processing is an option to 
consider for very large finite element problems. 

A comparison of analytical data with any experimental data will provide useful infor- 
mation into the layerwise analyses' strengths and deficiencies. For example, work by 
Boitnott, Johnson, and Starnes [144] included a nonlinear failure analysis of pressurized 
composite panels. The work in Ref. [144] compared experimental failure results with 
analytical failure results. The analytical analyses were conducted in order to simulate 
the actual experimental conditions. The analyses utilized the measured radius, 
circumferential slip, and axial strain for each experimental specimen to model the re- 
sponse as accurately as possible. The failure analyses of the curved panels described in 
Ref. [144] could also be accomplished using a layerwise finite element model, but it will 
involve a great deal of work and therefore will be left for future study. A study of 
geodesically stiffened shells subjected to high pressures may provide interesting 
interlaminar stresses. Additional analyses that incorporate geometric nonlinearity 
should be considered. Also, a study of a variety of laminate layups as a function of 
various shell and stiffener parameters may provide some unique results and a good da- 
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tabase. Additional failure theories and post-ply failure could easily be implemented and 
studied. The aforementioned recommendations coupled with the work accomplished in 
this research should provide several additional research projects. 
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Stiffness Terms 
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Appendix C 

Finite Element Stiffness Terms 
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C.l Layerwise Shell Element Direct Stiffness Terms 
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C.3 Layerwise Beam Element Direct Stiffness Terms 
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C.4 Layerwise Beam Element Tangent Stiffness Terms 
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C.5 Computation of Higher Order Derivatives 


The computation of the second and higher order derivatives of the interpolation func- 
tions with respect to the global coordinates involves additional computations. 


The first order derivatives with respect to the global coordinates are related to those with 
respect to the local (or element) coordinates according to 



dx 
d\l> n 

l“J 


dx_ 

H 

dx 

drj 



UT' 




(C.5. 1 ) 


where the Jacobian matrix [J] is evaluated using the approximation of the geometry: 


x - ]>\/^ (£,»/) 

J-i 

r 

y - 

j-i 


(C.5. 2) 


where <£/ are the interpolation functions used for the geometry and (f, rj) are the element 
natural coordinates. For the isoparametric formulation r = NDS and ft ■* ft. The 
second order derivatives of i ft with respect to the global coordinates (x, y) are given by 
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dy 2 
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^ dy dx 

34 <94 
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dy 54 d V 34 

A 

34 2 34 2 

d 2 x 3y 

3»f 2 d>/ 2 
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d s d i 
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d*l 34 


(C. 5.3) 


(C.5.4) 


(C.5.5) 


The matrices [ J\ ] and [ Ji ] are computed using Eq. (C.5.2). 
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